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Abstract 



\ The review is devoted to a relativistic formulation of the first Dirac quantization 

of QED (1927) and its generahzation to the non-Abehan theories (Yang-Mills and 
. QCD) with the topological degeneration of initial data. Using the Dirac variables 

\ we give a systematic description of relativistic nonlocal bound states in QED with a 

choice of the time axis of quantization along the eigenvectors of their total momentum 
operator. 

i We show that the Dirac variables of the non-Abelian fields are topologically de- 

generated, and there is a pure gauge Higgs effect in the sector of the zero winding 
number that leads to a nonperturbative physical vacuum in the form of the Wu-Yang 
monopole. Phases of the topological degeneration in the new perturbation theory are 
I determined by an equation of the Gribov ambiguity of the constraint-shell gauge de- 

fined as an integral of the Gauss equation with zero initial data. The constraint-shell 
non-Abelian dynamics includes zero mode of the Gauss law differential operator, and 
a rising potential of the instantaneous interaction, that rearranges the perturbation 
series and changes the assymptotic freedom formula. 

The Dirac variables in QCD with the topological degeneration of initial data de- 
scribe constituent gluon, and quark masses, the spontaneous chiral symmetry break- 
ing, color confinement in the form of quark-hadron duality as a consequence of sum- 
ming over the Gribov copies. A solution of U(l)-problem is given by mixing the zero 
mode with 770 - meson. We discuss reasons why all these physical effects disappear 
for arbitrary gauges of physical sources in the standard Faddeev-Popov integral. 

Is devoted to the memory of V.N. Gribov and I.V. Polubarinov 
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1 Historical Introduction 



The first papers by Dirac [jT|], Heisenberg, Pauli 0, and Fermi on the quantization of 
electrodynamics ran into difficulties of the determination of physical variables. The interpretation 
of all gauge components as independent variables contradicts the quantum principles; whereas 
excluding nonphysical variables contradicts the relativistic principles. 

The consistent quantum description of gauge constrained systems was considered by Dirac, 
Schwinger, Feynman, Faddeev, Gribov, and other physicists [jl|, ^, 0] as one of most 

fundamental problems of theoretical physics in the 20th century. 

The first quantization of electrodynamics belongs to Dirac who disregarded the relativis- 
tic principles and excluded nonphysical components by the reduction of the initial action to the 
solution of the Gauss law constraint, i.e. the equation for the time-like component of the gauge 
vector field. The Gauss law connects initial data of the time-like component with the data of 
all other fields. But this constraint- shell method had a set of defects, including nonlocality, 
explicit noncovariance as the dependence on the external time axis of quantization, and com- 
plexity. Feynman-Schwinger-Tomonaga formulation of QED admitted a simpler method based 
on the extended dynamics where all components were considered on equal footing with fixing a 
relativistic gauge. 

At the beginning of the 60s, Feynman found that the naive generalization of his method 
of construction of QED did not work for the non-Abelian theory. The unitary S-matrix in the 
non-Abelian theory was obtained in the form of the Faddeev-Popov (FP) path integral by 
the brilliant application of the theory of connections in vector bundle. There is an opinion that 
the FP path integral is the highest level of quantum description of gauge relativistic constrained 
systems. In any case, just this FP integral was the basis to prove renormalizability of the unified 
theory of electroweak interactions in papers by 't Hooft and Veltman awarded the Nobel prize 
in 1999. 

Nevertheless, in the context of the first Dirac quantization and its Hamiltonian generaliza- 
tions [^, |TU|, the intuitive status of the FP integral was so evident that two years after the 



paper (§[], Faddeev gave the foundation of the FP integral by the construction of the unitary S- 
matrix [^] for an "equivalent non-Abelian unconstrained system" derived by resolving constraints 
in terms of the radiation variables of the Hamiltonian description. 

Faddeev showed that on the one hand, the constraint-shell dynamics is compatible with the 
simplest quantization by the standard Feynman path integral, on the other hand, this Feynman 
integral is equivalent to the FP integral in an arbitrary gauge. This equivalence was proved by 
the change of variables in the Feynman path integral that removed the time - like vector of 
the canonical quantization into the phase factors of physical source terms. These phase factors 
disappear for S-matrix elements on the mass shell of elementary particles. In other words, 
Faddeev proved the equivalence of the constraint- shell approach with quantization of gauge 
theories by the gauge-fixing method only for scattering amplitudes where all color particle- 
like excitations of the fields are on their mass-shell. But the scattering amplitudes for color 
particles are nonobservable in QCD. The observables are hadrons as colorless bound states where 
elementary particles are off mass-shell. Just for this case, the Faddeev theorem of equivalence 
of different gauges becomes problematic even for QED in the sector of instantaneous bound 
states, as the FP integral in a relativistic gauge loses all propagators with analytic properties 
that lead to instantaneous bound states identified with observable atoms. 

The Faddeev equivalence theorem was proved before the revelation of nontrivial topo- 
logical properties [0, |13|] and the Gribov ambiguities 0. These new facts were studied 



mainly on the level of the FP integral. The topological degeneration of classical vacuum was 



incorporated into the FP integral in the form of the instanton solutions [jTlJ. 

It seems more natural first, to study the topological degeneration of all non-Abelian physical 
states at the more fundamental level of the constraint- shell dynamics compatible with the 
Feynman path integral and, then, to obtain the corresponding non-Abelian FP path integral 
analyzing possible Gribov copies. 

The first quantization of electrodynamics was fulfilled just at the level of the constraint- shell 
dynamics by Dirac |]T|. The Dirac reduction is based on the honest resolution of the Gauss law 
constraint and introduction of the Dirac nonlocal gauge-invariant physical variables |I4|, |l5l , 
instead of gauge-fixing. The Dirac method is the way to distinguishe the unique (radiation) 
gauge of physical sources. 

The present review is devoted to the generalization of the Dirac variables [0, |14|, ^ to the 
non-Abelian theories with the topological degeneration of the initial data in QCD in the class 



of functions of nontrivial topological transformations [|^, [T^, [T§, |1 

We present here a set of arguments in favor of that the constructed equivalent unconstrained 
system contains the most interesting physical effects of hadronization and confinement in QCD 
that can be hidden in explicit solutions of constraints and equations of motion P^, ITBI, 11^, ^21. 



We show also that the FP integral in an arbitrary relativistic gauge loses all these effects. 
The relativistic covariant properties of the Dirac variables allow us to use the Markov- 
Yukawa prescription of construction of the multilocal irreducible representations of the 
Poincare group and to formulate the generalized S-matrix formalism where the time axis of the 
Hamiltonian description is proportional to eigenvectors of the total momenta operator of any 



physical state [25, 26, 27 ||, 29, 3C, 31 



2 Dirac variables in QED 



2.1 Gauge- fixing method of 1967 

We formulate the statement of the problem using QED. It is given by the action 



W[A,i;,^] 



dx 



(1) 



d^-ieA^ )(7 = V^-7' 
d^A^ — djjA^ . 



(2) 



This action contains gauge fields more than independent degrees of freedom. First of all it 
is invariant with respect to gauge transformations 



A 



A^ + d^K ,7/^^ = exp[ieA]z^ . 



(3) 



One supposes that this invariance allows us to remove one field degree of freedom by the help 
of arbitrary gauges 

F{A^) = 0, FiA;) = MfUj^O, (4) 

where the second equation means that the gauge fixes field unambiguously. 

The standard quantization of the gauge field system in the gauge F(A^) = is based on 
the FP path integral 



l[DA^,Di^^Di;^^^FpS{F{A^))e 



(5) 



A 



where A^p = detMp is FP determinant and 

= j d^x [s^i;^ + i,^s^ + Ap^) (6) 

are the sources. The foundation of the intuitive FP integral by the canonical quantization was 
made in Faddeev's paper [|^. 

After Faddeev papers of 1967-69 there dominates very popular opinion that all physical 
results do not depend on a gauge F{A^) = 0. A gauge is defined for reasons of simplicity 
and convenience, including the condition detMp 7^ 0; the opposite case detMp = is called 
the Gribov ambiguity [^. All applications of quantum gauges theory after 1967-69, includ- 
ing investigation of the Gribov ambiguity |]32|], topological degeneration of initial data 



hadronization |Q, parton models [0], etc., were fulfilled at the level of the FP integral (H). 

While, in Faddeev's paper there was established the range of validity of the FP integral 
and its gauge independence. Faddeev managed to prove the FP integral by the canonical 
quantization only in the sector of scattering amplitudes for elementary particles on their mass- 
shell. This proof becomes doubtful for bound states where elementary particles are off mass-shell, 
the latter is crucial for QCD where only bound states are observable. 

The present review is devoted to consideration of the problems of bound states, Gribov 
copies, and topological degenerations of initial data at the more fundamental level of the Dirac 
canonical quantization of 1927 

2.2 Dirac reduction and Dirac variables 

The situation with canonical quantization of gauge theories before the FP revolution of 1967 was 
presented in the review by I.V. Polubarinov [|I4|. Igor Vasil'evich intended to sent his manuscript 
to Usp. Fiz. Nauk. It is pity that this nice review was not published. After 1967, Polubarinov 
tried to include the FP integral in his review, but he was not satisfied by the level of the physical 
foundation of the FP - scheme of quantization in the context of the first (constraint-shell) 
formulation of QED that belongs to Dirac [|l|. We reproduce here this formulation. 

The Dirac quantization was based on a definite frame of reference distinguished by the time 
axis T]^ = (1, 0, 0, 0), which allows us to establish boundary conditions and initial data. In special 
relativity, Einstein identified a frame of reference with a set of the watches, rulers, and other 
physical instruments for measurement of physical quantities (that include in our case a spectrum 
of the field excitations). The classical equations are split on the Gauss law constraint 

^ = ^ AAo = d^doAi + Jo (A = dA, = eV^T/.V') , (7) 

OAq 

and equations of motion 

6W „ 



5A, 



6W 

sw 



dlAk - dkdoAo - {5kiA - dkdi)Ai = , (8) 
= ^ ipii J7{A) + m°) = , (9) 



^ {i f{A) - m^)tlj = . 
oyj 

The problem of canonical quantization meets with the non-dynamic status of the time compo- 
nent Aq = A^T]^. The non-dynamic status of Aq is not compatible with quantization of this com- 
ponent as the fixation of Aq (by the Gauss law) and its zero momentum Eq = dC/d{doAo) = 



contradict to the commutation relation and uncertainty principle. To keep quantum principle, 
Dirac excluded the time component using the Gauss law constraint (|^): an explicit solution of 
this Gauss law 

Ao{t,x) = ao[A] + -jo{t,x) , (10) 

where 

ao[A] = -didoA{t,x) , (11) 

connects initial data of ^0(^072;) with a set of the initial data of the longitudinal component 
didoAi{t,y), and current joit,y) in the whole space, here 

= (12) 

A'' ^ ^ in J ^\x-y\ ^ ^ 

is the Coulomb kernel of the nonlocal distribution. One can substitute the solution (|T^) into 
equation for spatial component (^) 

6W 

[Sik - di—dk]idQ - A)Ak = ji - d,—dojo- (13) 



6A. 



i SAq 

We can see, that the constraint-shell equations ( |I3D contain only two transverse physical variables 
as the gauge - invariant functionals 

A*{t,x) = [6,k-d~dk]Ak . (14) 
Dirac generalized these gauge - invariant variables onto all fields by gauge transformations 

T.4A? = Am = v[A]{A, + z-d,){v[A])-\ 

a=l,2 ^ 

^^[AV^] = v[A]^, (15) 
where the gauge factor |]T|] is defined by 

v[A] = exp{ie dt'ao{t')} . (16) 

Jin 



Using the gauge transformations 

ao=ao + 9oA v[A^] = exp[ieA{to,x)]v[A]exp[ieA{t,x)] , (17) 

we can find that initial data of the gauge - invariant Dirac variables are degenerated with 
respect to stationary gauge transformations 

[A^] = [A] + dMto,x), ^""[A^,^^] = exp[zeA{to,x)]ij''[A^] . (18) 

The Dirac variables ( p3D as the functionals of initial fields satisfy the Gauss law constraint 

do{diA^{t,x))=0 . (19) 

Thus, explicit resolving the Gauss law allows us to remove two degrees of freedom and to reduce 
the gauge group into the subgroup of the stationary gauge transformations 



We can fix a stationary phase Aito,x) = ^o{x) by an additional constraint in the form of 
the time integral of the Gauss law constraint (^) with zero initial data 



= (20) 

We call this equation the constraint- shell "gauge". This "gauge" restricts initial data to a 
phase distinguished by the equations A$o(a;) = 0. Nontrivial solutions of this equation we call 
the degeneration of initial data, and the Gribov copies of the constraint-shell "gauge". The 
degeneration of initial data is determined by topological properties of the manifold of stationary 
gauge transformations in the class of functions with a finite density of energy. In the case of 
the three-dimensional QED there is only a trivial solution $o(a^) = 0. In this case 

ijj^ = ijj*, = A* . 

In the one-dimensional QED and non-Abelian theory the degeneration of initial data is the 
evidence of a zero mode of the Gauss law constraint that describes the topological excitation 
of gauge field with the Coleman spectrum of the electric tension [jT7|, |19|, We consider 

this zero mode in the Section devoted to Yang-Mills theory. 

Dirac constructed an equivalent unconstrained system the equation of which reproduce the 
equations of the initial theory (|l|). 



W* = WUw/SAo=o = J d'x^-[ E {d^Ald^AD + - J* A* + r{td - m)r] ■ (21) 



a=l,2 

To derive this equivalent unconstrained system that contains only physical variables, Dirac [Q, 



14| proposed to change the order of constraining and varying. He substituted the solution of 
the Gauss law constraint (0) into the initial action (|l|) in the rest frame r]^ = (1,0,0,0) and 
introduced the Dirac variables (p!5|). 

To combine the nonlocal physical variables A* and variational principle formulated for local 
fields, we can introduce three independent variables. In this case, the Dirac action should be 
added by the Lagrange multiplier ^ 

= W* + J d''xXL{x)diA* . (22) 

The local equations of the equivalent unconstrained system (|2^ ) 

SA* SXl 5ip* 5^* ^ ^ 

completely coincides with the equations of the initial constrained action (|^)- (^). Eqs. ( ^31) 
reproduce the constraint and lead to the equation for the Lagrange multiplier 

AAL(t,f) = 0. (24) 

The latter coincides with the equation for the stationary phase. In three-dimensional QED both 
they A, and $o si's equal to zero, in the class of functions where the physical fields are defined. 

In three-dimensional QED there are only three subtle differences of the equivalent uncon- 
strained system (^) from the initial gauge theory First of them is the origin of the current 
conservation law. In the initial constrained system (|lD, the current conservation law Sojo = diji 
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follows from the equations for the gauge fields; whereas the similar law doj^ = dij* in the equiv- 
alent unconstrained system ( ^if) follows only from the classical equations for the fermion fields. 
This difference becomes essential in quantum theory. In the second case, we cannot use the 
current conservation law, if the quantum fermions are off mass-shell, in particular, in an atom. 
What we observe in an atom? The bare fermions, or dressed ones (0)? Dirac supposed |^ that 
we can observe only gauge mvanant quantities of the type of the dressed i\e\ds. Really, we can 
convince that dressed f\e\ds (|TB|) as nonlocal functionals from initial gauge fields are invariant 
with respect to the time dependent gauge transformations of these initial fields (H). 

The gauge invariance with respect to the time dependent gauge transformations is the second 
difference of the nonlocal Dirac variables (|T^ from the initial fields of the constrained system (P 
with usual transformational properties with respect to gauge and Lorentz transformations. 

The gauge constraint diAi = 0, in the gauge-fixing method, is associated with the relativistic 
noncovariance. Whereas, the observable nonlocal variables ( |I5D depend on the time axis by the 



relativistic covariant manner. Polubarinov's review [|14|| was mainly devoted to the relativistic 
covariant formulation of the Dirac quantization |^. 

The gauge-fixing method and its terminology "the Coulomb gauge" do not reflect these 
three properties of the Dirac observables in the constraint- shell QED (pf) : the off-mass-shell 
nonconservation of the current, gauge invariance (^, and relativistic covariance. 

In fact, the term gauge ( ^0]) means a choice of nonlocal variables, or more exactly, a gauge 
of physical sources associated with these variables. 

This Dirac construction of a relativistic covariant equivalent unconstrained system can be 
generalized also on massive vector theories p7|. The generalization of the first Dirac quantization 



on the non-Abelian theory, including QCD [^, 0, ^ is the topic of the present review 



2.3 Relativistic covariance 

Relativistic transformations of the Dirac variables are discussed in details in the Polubarinov 
review |[T^ (see also 0). if we make usual relativistic transformations of the initial 



(25) 



fields Ai,AQ,ip with the parameter 

SlAk = ei(x-(9o' - XQdi')Ak{x') + e^Ao , 

Slip = ei{x'ido> - x'odi>)ij{x') + ^ek['yi,-fj]ip{x') , 

then the physical variables (|T5|) suffer the Heisenberg-Pauli transformations 

Al[A, + 6lA] - Al[A] = 6lAl + d,A , 

(26) 

^*[A + 6lA,^ + 6l^]-^*[A,^] = 6lij* + ieA{x')^* , (27) 

were 

A[A*,j*] = e,j^{doAl + dk^j*) . (28) 

These transformation were interpreted by Heisenberg and Pauli (with reference to the un- 
published note by von Neumann) as the transition from the Coulomb gauge with respect to the 
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time axis in the rest frame 77° = (1, 0, 0, 0) to the Coulomb gauge with respect to the time axis 
in the moving frame (see Fig.l) 



1,0,0,0) 



Al = 



lAl' = {v'-A) = 



1,0,0,0) 



rj' = r]o + 5lr]o 



Lorentz 
Frame 



Gauge 



New 

Lorentz 
Frame 



New 

Gauge 



Figure 1. 

These transformations correspond the "change of variables" 



(29) 



so that they became transverse with respect to the new time axis i]' (or, from the point of view 
of the "gauge-fixing" method of reduction, the transformations (^Sf) , (^61) ) correspond to the 
"change of gauge" . 

In result we got the relativistic covariant separation of the interaction on the Coulomb 
potential (instantaneous with respect to the time- axis r/^J and on the retardation. 
The Coulomb interaction has the covariant form 



W, 



c 



j d'xdS\j;{x)Vc{z^)j;{y)8{^ ■ z) 



(30) 



Here 



j* = eip*fiij* , = Zf,- ri^{z -r]) , z^ = {x- y)^ , 



(31) 



Voir) 



(32) 



Finite Lorentz transformations from the time axis r]^ ' to the time axis r]^ ' where constructed 
in paper []I4| using the gauge transformations 



(33) 



where t/2,1 = f (2)^(1), and f(2),t'(i) the Dirac gauge factors (|I6|) for the time axes t]^'^^ and t]^^^ 
respectively. 



Q 



2.4 Quantization and Feynman path integral 

The initial action (HI) is not compatible with quantum principles. The Dirac formulation of 
the equivalent unconstrained system keeps the quantum principles by the value of excluding 
the nonphysical components. We quantize the equivalent unconstrained system with gauge- 
invariant physical variables (^5]). The corresponding commutation relations 

z[doA:{^,t) , A*{y,t)] = {6,,-d.^d,)6'{^-y) , 

{r^M,r{y,t)} = 53(x-y) 

lead to the generating functional for Green's function of the obtained unconstrained system in 
the form of the Feynman path integral 

Z;[s,s*,J*] = Jl[DA*D^*Dij*e'^*'^'^''^*''^*^+'^*, (34) 
j 

with external source terms 

s- = /A(r# + ^.v + j-A-). (35) 

By the construction of the unconstrained system tW\s generating functional is gauge - invariant 
and relativistic covariant. Relativistic transformation properties of the quantum fields should 
repeat the ones of the Dirac variables (|T5|) as non-local functionals of the initial fields. As it 



was shown in papers |4[ |3^, ^ that quantum theory with the gauge - invariant 



Belinfante energy-momentum tensor on the constraint 

T^. = F^xF^ + ^j^[id, + eA,]^-g^,L + 



I 



+ -9,[^r;,V'], (36) 



4 

r^i. = ]p^iv\iv - gi,vi^ - glip. 

completely reproduced the symmetry properties of the "classical" theory (^5]), (P^l), ([ 

^e,[Mofc,^*] = + ieh\A\3M\ ^ofc = / d^Xx^T^^ - tTo^] . (37) 

This Lorentz transformation of the operator quantization means on the level of the Feynman 
path integral the change of the time axis 

Zl^{s\s\r\ = Z*[Ls\Lr,LJ*]. (38) 

This scheme of quantization depends on a choice of the time axis. If one choose a definite 
frame of reference with the initial time axis, any Lorentz transformation turns this time axis 
by the relativistic covariant manner. In this meaning, the constraint dynamics is relativistic 
covariant. Another problem is to find conditions when measurable physical quantities and results 
of theoretical calculations do not depend on the time axis (identified with a physical device). 
This independence exists only for scattering amplitudes of particles onto their mass-shell [0]. 



in 



In this case, one can say about the relativistic invariance of the scattering amplitudes of the 
local degrees of freedom. But, it is well-known that the Green functions (in particular, one- 
particle Green function) and instantaneous bound states depend on the time axis. In general 
case, measurable quantities in electrodynamics depend on the time axis and other parameters 
of a physical device including its size and energy resolution []T7|. 

If a nonlocal process depends on the time axis, one should to establish a principle of the 
choice of this time axis. 

We shall use the generalization of the Markov- Yukawa principle supposed in |jT^, ^ ^ 



4^ : in concrete calculations the time-axis is chosen to be parallel to a total momentum of any 
state. 

In particular, this choice and the nonlocal relativistic transformations ( ^Tj) remove all infrared 
divergences from the one particle Green function in the radiation variables [|T^, pll, 



i{27T)^5'^{p — q)G{p — q) = I d'^d'^y exp{ipx — iqy) < 0\Ttjj* {x)^p* {y)\0 > 



(39) 



G{p) = Goip) + GoipMp)Go{p) + Oia^ 
a f d'^q 



Goip) 



m\ 



Ti'^i J q^ + le 



5ij - Qi^Qjj liGoip + q)lj + 1oGo{p + qho^ 



a 

4:71 



n(p), 



where Il{p) is 



1 



m{3D + A)-D{^-m) + -{^-m) 



> -|- m) 

p2 



In 



1 1 
m — p 




2p2 



2p^ 



and is a ultra-violet divergence. The transition to another Lorentz frame p'^ = {p'q,p' ^ 0) is 
accompanied by the additional diagrams < 0\Tijj*(x)5Lijj*{y)\0 > induced by the transformation 
(^). As result in another frame we get the same relativistic covariant expression depending on 
the new momentum p' (see details in |42|). 



2.5 Gauge equivalence theorem and FP integral 

Thus, the constraint-shell generational functional ( P^ is relativistic covariant ([3^ ) and gauge- 
invariant by the construction. The main difference of this functional from the intuitive FP 
integral (H) is the information contained in the solution of the constraint, i.e., the electrostatic 
phenomena of the instantaneous interaction, including the Coulomb-like bound states. 

The Faddeev-Popov integral for the generating functional of Green functions in the gauge 
F(A) = can be obtained from the Feynman integral by two steps: *) a change of 
variables, and **) a change of sources. 

*) the change of variables is fulfilled by the Dirac factors (0), (0) 



v[A^]iA^,+t-^d,){v[A^]y 



(40) 



v[A^] =exp{ie-d^Af} 



(41) 
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This change introduces additional degrees of freedom and the FP determinant App of the 
transition to new variables of integration. These degrees are removed by additional constraints 
F{A) = 0. Finally, the constraint- shelUunct\ona\ Z* (^4]) takes the equivalent form of the FP 
path integral 



z*[s*,s*,r] 



(42) 



all electrostatic monopole physical phenomena that depend on the time axis are concentrated 
in the Dirac gauge factor v{A) that accompanies the physical sources s*, s*, J*. 



S* 



d^x ({v[A'^])-H*^^ + ij^{v[A^])-'s* + J*A*[A' 



(43) 



Thus, this change of variables corresponds to the rearrangement of the Feynman diagrams. In 
particular, after the change of variables the sum of the Coulomb kernel and transverse photon 
propagator 



/C^(J) 



7(1)^7(2) , 



converts into the identically equivalent sum of the Feynman gauge propagator /C^ and the 
longitudinal term /C^ : 

/C^(J) =/C^(J) + /C^(J) , 

where 



(44) 



/c-(j) = -[j«jf -j^v: 



tW 7(2)1 



go - ' 
7(2)- 



^'{J) = [(goJi'^)(goJi'^) - {q.J^iqjJr) 
**) The next step is the change of sources 

= j d'x(s^ij^ + i;^s^ + Alr 

In the result we get the original FP integral 



1 



q2(gg - q2) 



X{DAlD^lj^D^l,^A^pp5{F{A^))e 



(45) 



(46) 



without the Dirac factors. Changing the sources we lose the Dirac factor together with the class 
of the spurious diagrams that remembered the electrostatic phenomena and instantaneous 
bound states. One of these spurious diagrams is the longitudinal term /C^ in equation (0). 
This longitudinal term disappears only on the mass-shell (because of the current conservation 
law Jq ' go = Ji ' qi)- As we have seen before, the Dirac nonlocal gauge-invariant observable 
currents do not satisfy the conservation law for particles off the mass-shell, in particular, in 
bound states. 

Really, the FP perturbation theory in the relativistic gauge ( |46|) contains only photon propa- 
gators with the light-cone singularities forming the Wick-Cutkosky bound states with the spec- 
trum differing from the observed one which corresponds to the instantaneous Coulomb interac- 
tion The Wick-Cutkosky bound states have the problem of a tachyon and the probability 
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interpretation. These problems where solved by the quasipotential method that introduces 
instantaneous bound states, and the corresponding time axis as the Markov- Yukawa prescrip- 
tion ig. 



Thus, in QED, the fundamental constraint-shell functional ( |34l ) coincides with the FP inte- 
gral in the sector of scattering amplitudes for elementary particles on their mass-shell 
that is not actual for solution of the problem of hadronization in QCD. For description of the 
bound state sector in gauge theories, including QCD, we have three alternatives: 

1) the FP integral (|4|), 

2) its quasipotential approach to bound states at the level of the FP integral [^, and 

3) the constraint-shell functional ( ^41) added by the Markov- Yukawa prescription of the choice 
of the time axis. 

After Faddeev papers of 1967-69 solutions of all problems of gauge theories, including the 
description of bound states, the Gribov ambiguity 0, topological degeneration [H^, |I2|, were 
considered only at the level of the FP integral (|4^). 

The task of the present paper is to consider these problems at the most fundamental level of 
the Hamiltonian approach to quantization of gauge theories. It is based on a definite frame of 
reference, that includes a choice of the time axis, initial data, boundary conditions, normalization 
of wave functions, bound states, etc.. It is useful to recall the following words by Max Born 
about quantum theory ([^, p. 108): " The clue is the point that quantum mechanics does 
not describe a situation in an objective external world, but a definite experimental arrangement 
for observing a section of the external word. Without this idea even the formulation of a 
dynamical problem in quantum theory is impossible. But if it is acceptable, the fundamental 
indeterminacy in the physical predictions becomes natural as no experimental arrangement can 
ever be absolute precise." Following to Max Born one can say that the Hamiltonian description 
of any quantum system is determined by "a definite experimental arrangement for observing a 
section of the external word". If it is acceptable, the dependence of the quantum description of 
nonlocal processes on a frame of reference becomes natural as any experimental arrangement 
is included in this frame. The Hamiltonian method determines the energy spectrum of physical 
states. Therefore, a direction of the total momentum of any physical state is distinguished. In 
the case of bound states, the time axis is chosen along the total momentum of any bound state 



in the context of the Markov- Yukawa solution [|^, ^] of the problem of the relativistic 
covariance in QED. 



3 QED of bound states: spectrum and S-matrix 
3.1 Markov- Yukawa prescription 

One of the first definitions of the physical bound states in QED belongs to Lord Eddington: "A 
proton yesterday and electron today do not make an atom" It is clear that we can observe 
experimentally two particles as a bound state A4(x,y) at one and the same time 

M{x,y) = e''''''ij{z,)6{zo) , (47) 
where and are the total and relative coordinates 

^. = ^^^ , z, = ix-y), . (48) 



This principle of the simultaneity has more deep mathematical meaning p5], as the 
constraint of irreducible nonlocal representations of the Poincare group for arbitrary bilocal field 

Mix,y) = Miz\X) 

z,-^M{z\X) = , M{z\X)=M{x,y) . (49) 

This constraint is not connected with the dynamics of interaction and realized the Eddington 
simultaneity. 

The general solution of the irreducibility constraint ( |49D can be written in the form of the 
expansion of the bilocal field J^{z\X) with respect to "in" and "out" plane waves 



where a^-p^ are coefficients of the expansion. $-^^(2;^), '^-p^lz^) are the normalized amplitudes 
in the space of relative coordinates orthogonal to the total momentum of an atom Va 

(4). = - 'PA.i^) . (51) 
It is clear, that at the point of the existence of the bound state with the definite total momentum 



Vaim any instantaneous interaction (|30|) with the time- axis rj^ parallel to this momentum r]^ ~ 

Va^ 

V>.Miz\X) ~ VA,Miz\X) = ^^MiX\z) (52) 

is much greater then any "retardation" interaction [^]. It is just our principle of the choice of 
the time- axis of the Dirac quantization of a gauge theory. A time axis is chosen to be parallel 
to the total momentum of a considered state. In particular, for bound states this choice means 
that the coordinate of the potential coincides with the space of the relative coordinates of the 
bound state wave function in the accordance with the Markov- Yukawa prescription [|4| and the 
Eddington concept of simultaneity In this case, we get the relativistic covariant dispersion 



law and invariant mass spectrum. The relativistic generalization of the Coulomb potential is not 
only the change of the form of the potential, but also the change of a direction of its motion 
in four-dimensional space to lie along the total momentum of the bound state. The relativistic 
covariant unitary perturbation theory in terms of such the relativistic instantaneous bound states 
has been constructed in [Q. In this perturbation theory, each instantaneous bound state in 
QED has a proper equivalent unconstrained system of the Dirac quantization. The manifold 
of frames corresponds to the manifold of " equivalent unconstrained systems". In this case, the 
bilocal fields ( |60D automatically belong to the irreducible representation of the Poincare group 

By analogy, we introduce for the N- local field the total and relative coordinates 

1 ^ 

i=l 

which are connected by the identity 

N 



Then, the generalization of the Markov - Yukawa condition takes the form 



(54) 



Let Vfj, be the eigenvalue of the operator for the total 4-momentum, and r^^ be the unit vector 
in the direction V{r]^ ~ P^). Owing to the condition (|54l ) the N- local function 



being the Fourier transform of with respect to all coordinates, depends only on the 

transverse relative momenta 



N 



0. 



(55) 



1=1 



See also the generalization of the Markov- Yukawa condition for three-local and N- 
local [|0l cases. 



3.2 Effective Lagrangian of bilocal fields 

The constraint-shell QED allows us to construct the "bound state" relativistic covariant pertur- 
bation theory with respect to "retardation" [^, Our solution of the problem of relativistic 
invariance of the nonlocal objects is the choice of the time axis as a vector operator with eigen- 
values proportional to total momenta of bound states ( ^2]) |]25[| . In this case, the relativistic 
covariant unitary S-matrix can be defined as the Feynman path integral 



Z;[s,s\J*] =<*\ D^*D^*e''^c 



jw*[i^*,ip*]+is* 



* > 



where 



j 

is the averaging over transverse photons 

Wc[-^,ip] = f d^x[ilj{x){ip - ie/i* - m°)^(x) + 



(56) 
(57) 



+ IJ dSmy)^{x))lC^^\z^ I X)(V^(x)^(y))]. 



(5^ 



Here 

^ = 9^'7^, IC^'i) is the kernel 

lC^'^\z^ I X) =fiV{z^)5{z-r])fi, 



ifl = V'^lf. = 1-V, = ^1^- Vf.{z ■ T])) , (59) 



and z and X are the relative and total coordinates defined in equation (|4^ ), and V{z-^) is the 
potential depending on the transverse component of the relative coordinate with respect to the 
time axis t] 



d 



dX„ 



(60) 



In the context of this perturbation theory the bound state total momentum operators (23) 
can form the "new quantum numbers" of the type the Isgur- Wise ones |5^ . 



It seems that a most straightforward way for constructing a theory of bound states is the 
redefinition of action ( ^81) in terms of bilocal fields by means of the Legendre transformation 



= -l I d'xd^yMix,y)IC-\x,y)Mix,y)+ (61) 

+ / d^xd^y{ilj{x)ip{y))M{x,y) 

where JC^^ is the inverse of the kernel (0). Following ref. we introduce the short - hand 
notation 

d'^xd^yilj{y)ilj{x){ip - iej\* - m°)6^*\x - y) = {iPtfj , -G^^) , (62) 

d'^xd^y{i;{x)ip{y))Mix,y) = i'ipip,M). (63) 



After quantization ( or integration ) over Nc fermion fields and normal ordering, the func- 
tional (1^) takes the form 

Z?[s,s*,J*] =<*\ Jl[DMe'^^ff^^^+'^^ff^^^\*> (64) 

where 

W,jf[M] = (^^, i-G^' + A^)) - ^{M, IC-'M) (65) 
is the effective action, and 

S,ff[M] = {s*s*,{G-j,^-M)-^) (66) 
is the source term. The effective action can be decomposed in the form 

1 °° 1 

Weff[M] = --N,{M, IC-'M) + iN, (67) 

^ n=l 

Here $ = Ga-M.-, etc. mean the following expressions 

^{x,y) ^GaM = J d'^zGAix,z)M{z,y), 

= J d^xd'^y^{x,y)^{y,x), (68) 
$^ = y" d'^xd'^yd'^z^{x,y)^{y, z)^{z,x) ,etc 

As a result of such quantization, only the contributions with inner fermionic lines ( but no 
scattering and dissociation channel contribution ) are included in the effective action since we 
are interested only in the bound states. 

The requirement for the choice of the time axis ( |5UD in bilocal dynamics is equivalent to 
Markov - Yukawa condition ^ (||) 

'^'' ^^dxP = (69) 

where z^ = (x — y)^ and = (1/2) (a; + y)^ are relative and total coordinates. 



3.3 Quantization of bilocal fields 

The first step to the quantization of the effective action is the determination of its minimum 



n: 



5M 



n=l 



G-/-M 



0. (70) 



We denote the corresponding classical solution for the bilocal field by S(a; — y). It depends only 
on the difference x — y because of translation invariance of vacuum solutions. 

The next step is the expansion of the effective action around the point of minimum = 

^ + M' , 



W^fAM') = Wq{J:)+N,[-^M'IC-'M' + -iGj:M 



eff 

OO OO 1 

n=3 n=3 ^ 



iG^ = iG-/-J:)-% (71) 



and the representation of the small fluctuations Ai' as a sum over the complete set of orthonor- 
malized solutions F, of the classical equation 



6M^ 



\M'=0 



r = 0. 



(72) 



with a set of quantum numbers {H) including masses Mh = \/Pf and energies ujh 



M'{z\X)=Y. 



H 



{e^^^rH(g^|P)a+ (P) + e-^^^f H(g^| - V)a],{V)} 



The bound state creation and annihilation operators obey the commutation relations 



a%{V),a%,{V') 
The corresponding Green function takes form 



. 



(73) 



(74) 



(75) 



(76) 



^ {Vo - ujh - ie)2ujH {Vq - uh - ie)2ujH ' \ 
To normalize vertex functions we can use the "free" part of effective action (|7l| ) for the 



quantum bilocal meson M with the commutation relations (|74l). The substitution of the off - 
shell \/V^ 7^ Mh decomposition (|73| ) into the "free" part of effective action defines the reverse 
Green function of the bilocal field Q{Vq) 



W^%[M] = 27r5(Po -n)E / ^<i^) <^h{-'P)G-h\Vo) 



(77) 
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where Gh (^o) is the reverse Green function which can be represented as difference of two terms 



where M^{uj) is the eigenvalue of the equation for small fluctuations (^) and 



/(VP2) 



(27r) = 



(79) 



According to quantum field theory, the normalization condition is defined by formula 



2uj 



dVo l^°="(^^) dVo dM 
Finally, we get that solutions of equation (^) satisfy the normalization condition ^5 



dM{Vo) dI{M) 



\Vo=uj ■ 



iNr 



d 



dn 



dS 



tr 



and 



0^ (2vr)4 



2uj 



H 



^0) 
^2) 



is the fermion Green function. 



3.4 Schwinger-Dyson equation: the fermion spectrum 

The equation of stationarity ( [70D can be rewritten form of Schwinger-Dyson (SD) equation 

S(x -y)= m%^^\x -y) + ilC{x, y)G^{x - y) . (83) 
It describes the spectrum of Dirac particles in bound states. In the momentum space with 

S(A;) = j c/^xS(x)e*'=^' 

for the Coulomb type kernel we obtain the following equation for the mass operator ( S ) 

d\ 



4) 



(2 

where Gj:{q) = (^ — S(g))~^, l^(A;^) means the Fourier transform of the potential, 
kfj. — '>l^i{k ■ 7]) is the transverse with respect to r/^ relative momentum. The quantity S depends 
only on the transverse momentum 

S(fc) = S(fc^), 

because of the instantaneous form of the potential V_{k-^) in any frame. The fermion spectrum 



can be obtained by solving the Schwinger-Dyson equation ( |83|) . We may integrate it over the 
longitudinal momentum qo = {q ■ i]) using the representation 



^{q)=4^ + Ea{q^)S-\q^) 



i5) 
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for the self - energy with 

S-'{q^) = exp{-i^2va{q^)l = g^Vl^^l (86) 

where Sa is the Foldy - Wouthuysen type transformation matrix with the parameter Va- 
Then, one has 



where 



qo-Ea{q^)+ie qo + Ea{q^) + ie 



ASS.(g^) = 5a(g^)ASS(0)5-i(g^), a[1)(0) = (1 ± f,)/2 (88) 



]?) (87) 



are the operators separating the states with positive ( +Ea ) and negative ( —Ea ) energies. As 
a result, we obtain the following equations for the one - particles energy E and the angle v : 

E,{k^) cos2v{k^) =K + \j ^H^^ - q^) cos2v{q^) (89) 



Ea{k^) sm2v{k^) = \k^\ + - ■ sin2v{q^) (90) 

3.5 Bethe - Salpeter equation: bound-state spectrum 



Equations for the spectrum of the bound states ( [72D can be rewritten in the form of the Bethe- 
Salpeter (BS) one 



r = i}C{x, y) j d!^ zid^ Z2G-^{x - Zi)T{zi, Z2)G-^{z2 - y) . (91) 
In the momentum space we obtain with 

T{q\V) = J d^xdSe'^^d^^~y'>''T{x,y) 
for the Coulomb type kernel the following equation for and the vertex function ( £ ) 

where V_{k-^) means the Fourier transform of the potential, /c^ = k^ — r]^{k-r]) is the transversal 
with respect to r]^ relative momentum, is the total momentum. 
The quantity £ depends only on the transversal momentum 

mv) = T{k^\v), 

because of the instantaneous form of the potential V_{k^) in any frame. 



1Q 



We consider the Bethe - Salpeter equation ( |91|) after integration over the longitudinal mo- 
mentum go- The vertex function takes the form 

Ta,{k^\V) = J ^nk^ - q^Wa,{q^)l (93) 

where the bound state wave function 'ipab is given by 

^ ^^ A(+)a(9^r,,(g^|P)A(_),(g^) ^ A(-)a(g^r,,(g^|P)A(+),(g^) ^^ ^^^^ 
Et — VV^ + ie Et + VV^ — ie 

Et = Ea + Eb means the sum of one - particles energies of the two particles (a) and { b ) 
defined by (|9D, (|0D and the notation 

A(±)(g^) = 5-i(g^)A(±)(0)5(g^) = A(±)(-g^) (95) 

has been introduced. 



Acting with the operators (|88| ) and ( |95D on equation one gets the equations for the 



wave function ip in an arbitrary moving reference frame 



Erik^) T VV^)A'(i^Sk^)^Pab{k^)A'l^l)bi-k^) = 

^^&aik^)[J ^,nk^ - <l^)^ab{q^)]^^;^^.{-k^). (96) 



All these equations ( P3D and (^) have been derived without any assumption about the 
smallness of the relative momentum |A;^| and for an arbitrary total momentum 



We expand the function \& on the projection operators 

^ = ^+ + ^_, ^± = a|^^^A^^ (97) 
According to Eq. (p4D, satisfies the identities 

^ ^07) = , (98) 

which permit the determination of an unambiguous expansion of \& in terms of the Lorentz 
structures: 

^±(a,6) = 5(-4 {75l^±(a,6)(g^) + (7m - //)iVia,6)} A?(0)5(-,J , (99) 

where L± = Li±. L2, N± = Ni ± N2. In the rest frame r^^ = (1, 0, 0, 0) we get 

N>^ = (0, N') ; N\q) = ^ N^{q)eM + S(g)g^ . 

a=l,2 

The wave functions L, N", E satisfy equations 

1. Pseudoscalar particles. 

90 



r da o 



M 



(100) 

Li = EL2+/^V^(p-q)(c+c+ + es+s+)L2 • 
2. Vector particles. 



/(id ^ 

^^(P - q){(s c-r'' + s^-s-(r''e - r^V)) iVi +(r]"s c+) Si} ; 



(101) 



77° = gA"(p) , t=PCet{q) , = er(g)ef (p) 
3. Scalar particles. 





MS2 = + 







(102) 



M El = + 

r da 0^0^ 

+ J (^^(p - q){(^s s + ^;^.") ^2 +(!z^c+c-) iV2 } 

Here, in all equations, 

c^{p) = cos[va{p) ± Vb{p)] , s^{p) = sm[va{p) ± ^^^(p)] (103) 

Va,viy are the Foldy-Wouthuysen matrices of particles (a,b), and E = E"- + \s the sum of 
one-particle energies. 

The normalization of these solutions is uniquely determined by equation 

^^'^ E (?) + LMLt{q)] = 1 , (104) 



Ml ^ (2 



^ E (0 {ivr(?)iV2"'"(?) + iV2"(?)ivr (?)} = 1 , (105) 

^ E {Si(g)s^(g) + S2(g)s+(g)} = 1 . (106) 

The description of instantaneous relativistic bound states in a hot and dense medium can be 
found in 0, g. 
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3.6 Schrodinger equation 

If the atom is at rest ( = (M^, 0,0,0) ) equation (^) coincides with the Salpeter equation 
|]55[| . If one assumes that the current mass mP is much larger than the relative momentum 
then the coupled equations (^) and ( |96D turn into the Schrodinger equation. In the rest 
frame ( Vq = Ma ) equations ( |8^ ) and (^) for a large mass (m°/|g^| oo) describe a 
nonrelativistic particle 

Ea{k) = v/K)2 + k2 ^ mO + 

tan2i; = ^ 0; S{k) ~ 1; A(±) ~ 

Then, in equation only the state with positive energy remains 

ip ~ ^(+) = A(+) 75^4/1^^5^/1, A(-)ipA(+) ~ 0, 
where /i = ■ mh/{ma + mi,). And finally the Schrodinger equation results in 

[_Lk-2 + (^0 + ^0 _ M^)]V^5ch(k) = / ^Z(k - q)^5ch(q), (107) 

with the normalization / d?q\il)sch\^ / {^iiY = 1- 

For an arbitrary total momentum equation ( |107|) takes the form 

f"^^^'^^^' + ^""^ + ~ v^)]V^5c.(A:^) = / ^,nk^ - Q^)^ScHiq^), (108) 

and describes a relativistic atom with nonrelativistic relative momentum \k^\ << In the 

framework of such a derivation of the Schrodinger equation it is sufficient to define the total 
coordinate as X = (x + y)/2, independently of the magnitude of the masses of the two particles 
forming an atom. 

In particular, the Coulomb interaction leads to a positronium at rest with the bilocal wave 
function 

= Vpit) {^l^y tsJ^)^/f , tsJ^) = J (|^e(--^^,,.(p-); (109) 
where ip^^f^iz) is the Schrodinger normalizable wave function of the relative motion 

^ ^ ^]^.J^) = ^^..M^ ( [d'z\\^,Az)f = l) (110) 



\ TTlf, dz^ \^\ J —Sch^ \ J —Sch 

Mp = (2me — e) is the mass of a positronium, is the projection operator on the state with 
positive energies of an electron and positron. We have chosen the total-motion variable r]p{t) 
in (|109| ) so that the effective action for the total motion of the positronium with anomaly term 



have the form similar to the tiq -meson ones [§9 



Weff = Jdt[l{vl- Mlr^p^) V(3) + Cpr^pX[A]^ , (111) 



9.9. 



where 



and 



Cp = —8^ \ , (112) 



-X[A] = J £xF^,*F'^^ ' X[A] = ^J d'x{e,,kAd,Ak) (113) 

is the "winding number" functional (i.e., anomalous term) that describes the two 7 decay of a 
positronium. 

3.7 Spontaneous chiral symmetry breaking 

The solution of the set of SD and Salpeter equations (|89| ), ( |90D , ( |96D was considered in the 
numerous papers 0, ^ (see also review ||5^) for different potentials. One of the main 
results of these paper was the pure quantum effect of spontaneous chiral symmetry breaking for 
non-Coulomb potentials. The instantaneous interaction, in this case, leads to rearrangement of 
the perturbation series and strongly changes the spectrum of elementary excitations and bound 
states in contrast to the naive perturbation theory. 

To demonstrate this effect and estimate possibility of the considered relativistic equations, 
we consider the opposite case of massless particles, m° = m° for an arbitrary potential. 
Suppose that in this case equations (|89| ), ( |90|) 

2Ea{k^) cos 2v{k^) = J ^^^V{k^ - q^) cos 2v{q^) (114) 

2Ea{k^) sm2v{k^) = \k^\ + [ ^^V{k^ - q^)\k^ ■ sm2v{q^) (115) 

J (zvr)'^ 

have a nontrivial solution v{k-^) 7^ 0. This solution describes the spontaneous breakdown of 
chiral symmetry [||, |31|, |5|, ^, |] 



It can easily be seen that equations ( p.l4| ) and ( |115D are identical with {J^^ for the bound 
state wave function with zero eigenvalue, = and 

A(+)^A(_) = A(_)^A(+) = ip 

2E,{k^)^{k^) = j ^V{k^-q^)i,[q^). (116) 

Therefore, 

il) = cos2v{k^) / F (117) 
where F is a proportionality constant determined by the normalization ( p.04| ) 

^ AN^ f d^q 



^/(d^^^^°^^2.(g)) (118) 



MhJ (27r) 

In this way, the coupled equations (pQ]), (0), and (^) contain the Goldstone mode that ac- 
companies spontaneous breakdown of chiral symmetry. Thus, in the framework of instantaneous 
action we get the proof of the Goldstone theorem in the bilocal variant. 



Just this example represents a model for the construction of a low - energy theory of light 
mesons, in which the pion is considered in two different ways, as a quark - antiquark bound state 
and as a Goldstone particle. So it turns out that our relativistic instantaneous model for bound 
states can, in the lowest order of radiative corrections, also describes mesons. It was shown that 
the spontaneous symmetry breaking is absent for the pure Coulomb potential in QED The 



spontaneous symmetry breaking and the Goldstone meson in QCD are realized for the potential 
of hadronization. What is an origin of this potential in QCD? 

Recall that in QED the Coulomb potential is the consequence of the resolving the Gauss law 
constraint. This Coulomb potential is lost in any relativistic invariant gauge of sources in FP 
representation of the path integral. The defect of the pure relativistic artificial gauges is problem 
of a tachion and other non-physical states in the bound state spectrum [^. The incorporation 



of the simultaneity for these gauges by the quasipotential approach could not describe both 
the spontaneous chiral symmetry breaking with the Goldstone meson and S-matrix elements of 
interactions of some instantaneous bound states (papers on this topic are absent in literature). 

3.8 The relativistic equation for multiparticle systems. 

To derive the relativistic covariant equations for many - particle, we use the operator approach 
5D| , ^ with the Hamiltonian given by 



n = J d^^itdar + m°)^ + i I rfxrfy(v^+(x)v^,(x))r(x - y)(v^,+(y)^Ky)) • im 

The first step for constructing the physical states consists in the definition of the one - quasi 
- particle creation ( a+,6+ ) and annihilation {a,b) operators with the help of the Bogoliubov 
fermion expansion 

M^) = E / (^^e^^"K(q)/^-(q, s) + bt{-q)u4-q, s)] . (120) 

Here /iQ,(q, s) and z/q,(— q, s) are the coefficients determined from the Schrodinger equation for 
the one - particles energy 

< a,{q)\H\at{q') >= E{q) < 0|a,(q)a+(q')|0 > . (121) 

They can be represented via the Foldy - Wouthuysen matrix ( P5| ) as 

/i«(q, s) = S{q)af3fxp{0, s); //^(-q, s) = S{-q)afii^(3i0, s) 



with 



5„.Kq)[E/^«'(0,3)/.+ (0,.)]5,-,;,(q) = {S^-^S~') ^ (A°(q))^ 

s Z a/3 



/3' 



Saa'i-q)[j:^Ao,s)4io,s)]s,,li-q) = is^--^s-') ^ (A°(-q)).,- 

A° and A° are projection operators on states with positive, resp., negative energy. Then, 
equation ( |121D takes the form of the Schwinger - Dyson equation (|5^), (|9D| ) which can compactly 
be written as 

Eip)S-\p) = m'+pa^ + \j (I^Hp - ^)S-\q) , (122) 
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After inserting ( |120|) into ( |119|) the Hamiltonian can be given in the following manner 



7i = Eq -\- Hi+ : H4 :, 
Eo = <0|7^|0>, 

H, = ^E(pi)(a+ai + 6+6i), (123) 

(1) 
2 

■ H^: = - ^ 5^^) {pi -P2+P3- P4) V:(pi - Ps) 

1,2,3,4, 

+ 6ia2aj]"6| z/|/i2/i3^'4 + &ia263a4Z/j/i2Z^|/i4 + ...} + ... 
The following abbreviations have been used in ( |123D : 

E = E / {1} = {Pr^ ^/}. = {-P/, -^/}, / = 1, 2, 3, 4. 



For diagonalizing the Hamiltonian ( p.23| ) with respect to pair correlations {afb^), {b-^a^} one 
defines a new vacuum as the coherent state 

|0 »,= exp{ Y: «(l,2,3,4)[(a+^^6+^^)(6+^'^a+^'^)]}|0 > (124) 
1,2,3,4 

and the creation operator for the bound state (of pair correlation) 

S+(n)=^5(pi-p2)[X+(l,2)a+Xl)&+^(2)-X_(i,2)6^(i)a^(l)]. (125) 
1,2 

The coefficient X+ and X_ are determined from the Schrodinger equation for the two - particle 
energy Mb , 

„ << 0\B{n){Hi + Hi)B+{n)\0 »„= Mb a « 0\B{n)B^{n)\0 »„, (126) 



and the parameter a in ( |126|) is given with the help of the definition of the annihilation operator 



B{n) for the pair correlation 

B{n)\0»a=0. (127) 

Equation ( |126|) coincides with equation (|96|) in the rest frame (the Salpeter equation) for the 
meson spectrum 

{E,{p) + E^ip) T Mb)^±±{p) = A±(p)[/pq X (7/>++(q) + V^__(q))]A±(-p) (128) 
up to the notation 

ij = + i>—; V'±± = Aiz/'A:^; 

V'++(p)a/3 = E ^+(P' P' ^1' ^2)/ia (P, Sl)z//3(p, S2), (129) 

Sl,S2 

^— (p)a/3 = E ^-(P' P' ^1' ^2)z^a (P, Sl)flf3{p, S2), 



4q X Hq) = I (0i^(p - q)^(q) 
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The one - particle energies -E'i(p), E2{p) in ( p.29| ) are defined via the Schwinger - Dyson 
equation ( |122|) . 

Notice that equations of the type ( |122|) , (|129|) are well - known from the nonrelativistic 
many - body theory (Landau's theory of fermi liquids |0, Random Phase Approximation |]65[| ) 
and play an essential role in the description of elementary excitation in atomic nuclei [^. Their 
relativistic analogies describing the Goldstone pion and the constituent masses of the light quarks 
are equations O, (^. 

Thus the Green function method discussed in sects. 3. 2.-3. 5., and the operator approach lead 
to one and the same equations and complement each other. The first allows one to make easily 
the relativistic generalization and to construct the effective bound state interaction Lagrangian, 
whereas the second yields an adequate interpretation of quantum states and enables one to 
describe more complicated system, (in QCD, baryons and other many - quark states Pm ). 

Let us construct by means of the quasiparticle operator method the relativistic equation for 
a three - particle system. In the "coherent" vacuum ( p.25| ) the creation operator of a three - 
particle system consists not only of creation operators for particles (a+) but also of annihilation 
operators for antiparticles (6) with the same quantum numbers 



= E '^(Pi + P2 + P3)[X+++(1, 2, 3)a^(+)(l)a^-(+)(2)a'(+)(3) + 

1,2,3 

+ X^^+{l,2,3)b'^+\l)b''^+\2)b''^+\3)+ interchange of (l,2,3)]e*^^ (130) 
The "baryon" functions are as follows: 

V^+++(l,2,3)„^^ = E /^^(l)/iM2)/i+(3)X+++(l,2,3), 



S1S2S3 



^—+(1, 2,3)^^3^ 



E z.+ (l)z.+ (2)z.+ (3)X__+(l,2,3), 



S1S2S3 



etc. 

Then, the eigenvalue equation for the Hamiltonian operator 

„ «0\BHB+\0 

is equivalent to the following system for the "baryons" wave functions 



(131) 



E(2) 



+ 



3 



- / 


V 




V 


+ / 


1 + ] 


(1)A 


1 + ] 


(2)A 


+ 


V + ) 




+ 

I - J 




[ - 1 



+ 
+ 
V + / 



+ + + X (1,2,3) 



V +-- / 



(3) 



{ W^/ + + + N (1,2,3) . (1,2,3)] + 

+ + + 



V +-- / 



V - + - / 
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+ 43[^ / + + + . (1,2,3) + ^ . . (1,2,3)] + 



+ 



V +-- / 



+ - 



V + + + / 



+ hAi' I + + + s (1,2,3)+^, ^_ . (1,2,3)]} 



V +-- / 



+ 



V --+ / 



where 



7^^Hq)^(pi - q, P2 + q, ps) 



(2 

Pl + P2 + P3 = 0. 



(132) 



(133) 



Equation ( |133|) is the analogue of the Salpeter equation ( p.28| ) for a bound state consisting of 
three particles. In the same notations eq. ( [I29| ) has the form 



Mb]^±±(1,2) 



= -A±(l){42[^++(1,2) +^_.(1,2)]}A±(2) 
with the condition Pi = P2 = P and with the taking into account the identities 

/ dqV{p - q)V^(q) = 
= / dqV{q)ij{q + p) = 

= I c?qK(q)V'(pi + q, -P2 - q) \p,=p2=p 

^ (Pi,-P2)=^(l,2) 



The nonrelativistic reduction [^^l ^fom the Salpeter eq. ( |96| ) to the Schrodinger equation 



Eaip) 



T, — 2 , 1 P^ 

2 rrin 



Sa{p) ^ 1, ^+++ = ^»^ 



I + 



V 



leads in our case to the well - known nonrelativistic equation for the wave function of three 
particle bound states 



2 

o[4,2V^(Pl'P2'P3) +%3^(Pl'P2'P3) +43^(Pl'P2>P3)]- 



(134) 



Here, the condition (|134|) , which means the choice of the rest frame = (M^, 0,0,0) , has 
to be fulfilled. 
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Notice that the Jacobi coordinates, which allow to write the Hamiltonian in the term of two 
relative momenta, have sense only in the nonrelativistic limit. 

To describe the three-particle bound system in an arbitrary reference frame it is sufficient 
to substitute in (|132| ) all relative momenta by the transversal ones, p^*-*-*, and the projection 
operators A-|-(p) by the operators 



In the same way one can generalize the equation (|132|) and its relativization for an arbitrary A^- 
particle state. 

The method for constructing relativistic wave functions of many - quark system explained 
above unambiguously enables one to build from the nonrelativistic bound state wave function 

Aai,a2, ■■■ «jv ^ai,02, ■■■ Oat Vt^ 5 ••• /■> / , t-'V^/ " 

i 

relativistic wave functions for the same bound states with total momentum 



Here Xai,a2, ... ajv the matrix selecting one or another representation of the Lorentz group 
with a definite spin. ( A representation of the Poincare group that preserves the one - time 
dependence of wave functions see in ref. [|67| ). 

3.9 Relativistic covariant unitary S-matrix for bound states 

The achievement of the relativistic covariant constraint-shell quantization of gauge theories is 
the description of both the spectrum of bound states and their S-matrix elements. 

It is convenient to write the matrix elements for the action (0), ( [71] ) in terms of the field 
operator 

^\x, y) = J d'^xiGj:{x - xi)M'{xi,y) = ^'{z\X) 
Using the decomposition over the bound state quantum numbers (H) 



where 



H 



(135) 



^H(ab){q^\V) = G^a{q + V/2)TH(ab)iq^\V) , (136) 

we can write the matrix elements for the interaction (^) between the vacuum and the 

n-bound state 

Cn \ " 

E^O n 
1=1 ) j=i 



(27r)32cjj 



1/2 



M(")(Pi,...,P„) 

(137) 
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^ ^ {h} 



«3 



{{ik} denotes permutations over ik). 

Expressions ([73]), ([76|), ( |135|) , ( p.37| ) represents Feynman rules for the construction of a 
quantum field theory with the action (0) in terms of bilocal fields. 

It was shown that the separable approximation of the constraint-shell gauge theory of 



bound states leads to the well-known Nambu Jona-Lasinio model Q and the phenomeno 



logical chiral Lagrangians ^ used for the description of the low-energy meson physics. 
Thus, the constraint-shell gauge theory of bound states is sufficient for describing the spectrum 
and interaction of hadrons as extended objects (without introducing the ideology of bags and 
string). 

In the context of the constraint-shell gauge theory, to solve the problem of hadronization in 
QCD, one needs to answer the questions: 

i) What is the origin of the potential of hadronization in the non-Abelian theory? 

ii) How to combine the Schrodinger equation for heavy quarkonia (that is derived by the residuum 
of poles of the quark Green functions) with the quark confinement [[71] ? 

iii) What is the origin of the additional mass of the ninth pseudoscalar meson [|72| ? 



:i38) 



4 Dirac variables in Yang-Mills theory with the topo- 
logical degeneration of physical states 

4.1 Constraint- shell radiation variables in perturbation theory 

We consider the Yang-Mills theory with the local SU(2) group in four - dimensional Minkowskian 
space - time 

W[A,] = ~y d'^xF^^^Fr = \jd'x (Fof - Bf 
where the standard definitions of non-Abelian electric tension Fq^ 

Fo, = doA1-D{A)fAl Df = (5'^% + ge'^''^ A\ 

and magnetic one 

5f = e,,, {d,Al + ^-e^''^A]Ai) 
are used. The action ( p.38| ) is invariant with respect to gauge transformations u{t,x) 

A1 ■=u{t,x)(^Ai + di)u-\t,x), i)"" ■.= u{t,x)i) , (139) 



where A^ = g^A"^ . 

Solutions of the non-Abelian equations 



6W 







D^A) A', = D''\{A)doA'r (140) 
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5,,L)i {A) - D, {A) A {A) A] = o [A) d^A^ - D{A)fA', (141) 



are determined by boundary conditions and initial data. 

The first Gauss equation ( |140|) is the constraint. It connects initial data of Aq with the 
ones of the spatial components A°;. To remove nonphysical variables, we can honestly solve this 
constraint in the form of the naive perturbation series 



A'o = <^o[A] = ^dod,A1 + ... 



(142) 



The resolving of the constraint and the substitution of this solution into the equations of motion 
distinguishes the gauge - invariant nonlocal (radiation) variables. After the substitution of this 
solution into ( |141|) the lowest order of the equation (|141|) in the coupling constant contains 
only transverse fields 



(^0 - mi 



0. Af = [6,,-d^-dk]Al + ... 



(143) 



This perturbation theory is well-known as the radiation (or Coulomb 0) "gauge" with the 
generational functional of the Green functions in the form of the Feynman integral 



aTi 



J Jllld^A'^'^d^E'^^] 



c=l 



X 



(144) 



with the constraint-shell action ( [I38| ) 

M/Jo) [A^, E"^] = W' 



swi 



(145) 



given in the first order formalism 

W' = JdtJ d'x{F^,E^ - ^[E^E^ + B^B^]} 



The constraint 



Df{A)Ef = 



is solved in terms of the radiation variables 



E^ = Ef'^ + d,a, = 0. 



where the functions a"" take the form 



(146) 



(147) 



(148) 



(149) 



The operator quantization of the Yang-Mills theory in terms of the radiation variables belongs 
to Schwinger [@] who proved the relativistic covariance of the radiation variables (|143|) . This 
means that the radiation fields are transformed as the nonlocal functional 



Al[A] = v^[A]{A, + d,)iv^[A]) 



-1 



Ai 



Al^T" 
2i 



(150) 



where the matrix v'^[A] is found from the condition of transversallity dkAj^ = 0. At the level of 
the Feynman integral, as we have seen in QED, the relativistic covariance means the relativistic 
transformation of sources ^1 . 



The definition ( |150|) can be treated as the transition to the new variables that allows us to 
rewrite the Feynman integral in the form of the FP integral p, | 
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Z^[/W,J'^^] = I Jl[[d'A^]6{d,A'r)Det[D,{A)d,] 



X exp |ziy[y4] - i J 



dHJf-Af[A]) 



(151) 



It was proved [|, |73|, ^ that on the mass-shell of the radiation fields the scattering amplitudes 
does not depends on the factor It remained only to reply question: Why does nobody 

observe these scattering amplitudes of mass-shell non-Abelian radiation fields? There are some 
possible answers on this question: the infrared unstability of the naive perturbation theory |^ 



[75] , the Gribov ambiguity, or zero of the FP determinant the topological degeneration of 
the physical states [jT7|, |19|, In any case, in the non-Abelian theory there are not observable 
physical processes for which the gauge-equivalence theorem is valid. 



4.2 Topological degeneration of initial data 

One can find a lot of solutions of equations of classical electrodynamics. The nature chooses two 
types of functions: the monopole that determines nonlocal electrostatic phenomena (including 
instantaneous bound states), and multipoles that determine the spatial components of gauge 
fields with a nonzero magnetic tension. 

Spatial components of the non-Abelian fields considered above as radiation variables ( |143|) 
in the naive perturbation theory are also defined as multipoles. In the non-Abelian theory, there 
is the reason to count that the spatial components of the non-Abelian fields belong to the 
monopole class of functions like the time component of the Abelian fields. 

This fact was revealed by the authors of instantons [0]. Instantons satisfy the duality 
equation in the Euclidean space, so that the instanton action coincides with the Chern-Simons 
functional (Pontryagin index) 



u[A] 



tout 



167r2 



dt 



X[A^t] - X[Ain] = n{l 



out J 



n(ti 



(152) 



where 



X[A] 



d^xe'^^Tr 



AidjAk 



V 



—AiAjAk 



Ain,out — A{tin^Qut^X^ 



is the topological winding number functional of the gauge fields and n is a value of this functional 
for a classical vacuum 

i- = = . (154) 

The manifold of all classical vacua ( |154|) in the non-Abelian theory represents the group of 
three-dimensional paths lying on the three-dimensional space of the S'f/c(2)-manifold with the 
homotopy group 7r(3)(S'f/c(2)) = Z . The whole group of stationary matrices is split into topo- 
logical classes marked by integer numbers (the degree of the map) defined by the expression 

H\n\ = I £x e^^' Tr[L^L-^Ll] = n (155) 

which counts how many times a three-dimensional path v{x) turns around the S't/(2)-manifold 
when the coordinate Xi runs over the space where it is defined. 

Gribov in 1976 suggested to treat instantons as Euclidean solutions interpolating between 
classical vacua with different degrees of map. 

The degree of a map ( |155|) can be considered as the condition for normalization that deter- 
mines a class of functions where the classical vacua ( |154|) are given. In particular, to obtain 



Eq. ( p.55| ) we should choose a classical vacuum in the form 

t)(")(f) = exp(n<l>o(^)), $0 = -ivr /o(r) (r = |f|) , (156) 

r 

where the function foir) satisfies the boundary conditions 

/o(0) = 0, /o(oo) = 1 . (157) 

The normalization ( p.55| ) points out that the vacuum values L" of spatial components Ai belong 
to monopole-type class of functions. To show that these classical values are not sufficient to 
describe a physical vacuum in the non-Abelian theory, we consider a quantum instanton, i.e. 
the corresponding zero vacuum solution of the Schrodinger equation 

H^o[A] =0 {H= [ d^x[E'' + B\ E = ^) . (158) 

J 10 A 

It can be constructed using the winding number functional ( |153| ) and its derivative 

^XW = (159) 



The vacuum wave functional in terms of the winding number ( |153|) takes the form of a plane 
wave [ll6| 

^o[A]=eM^PNX[A]) (160) 



for nonphysical values of the topological momentum = ±i8n'^/g'^ |]T^, We would 

like to drew attention of a reader to that in QED this type of the wave functional belongs to 
nonphysical part of a spectrum like the wave function of an oscillator (p^ + g^)V^o = 0. The value 
of this nonphysical plane wave functional for classical vacuum (|154| ) coincides with quasiclassical 
instanton wave function 

exp(iVr[Anstanton]) = ^o[^ = Lout] X ^o[^ = ^m] = exp(-^[nout - ^in]) • (161) 



This exact relation between a classical instanton and its quantum version ( |158|) points out 
that classical instantons are also nonphysical solutions, they permanently tunneling in Euclidean 
space-time between classical vacua with the zero energy that does not belong to physical spec- 
trum 0. 

4.3 Physical vacuum and the gauge Higgs effect 



The next step is the assertion about the topological degeneration of initial data of not only 



classical vacuum but all physical fields with respect to stationary gauge transformations 

if) (to, x) = (to, x)t;(")(f)"' + , = v^''\x)div'^"\xy' . (162) 



The stationary transformations t>"(x) with n = are called the small ones; and those with 
n 7^ the large ones [^. 



The group of transformations ( |162| ) means that the spatial components of the non-Ahelian 
fields with a nonzero magnetic tension B{A) ^ belong to the monopole class of functions 
like the time component of the Abelian fields. In this case, non-Abelian fields with a nonzero 
magnetic tension contain a nonperturbative monopole-type term, and spatial components can 
be decomposed in a form of a sum of the monopole ^f' (x) and multipoles Ai 

Af\to,x) = ^f\x) + Af\to,x) . (163) 
The multipole is considered as a weak perturbative part with asymptotics at the spatial infinity 

Mto, ^) I asymptotics = C'(^) (/ > 1) • (164) 



Nielsen and Olesen [^, and Matinyan and Savidy ||7^ introduced a vacuum magnetic 
tension, using the fact that all asymptotically free theories are unstable, and the perturbation 
vacuum is not the lowest stable state. 

The extension of the topological classification of classical vacua to all initial data of spatial 
components helps us to choose a vacuum monopole with a zero value of the winding number 
functional (|153|) 

X[A = $f )] = , ^U=*(o) ^ 0- (165) 



The zero value of the winding number, transversality, and spherical symmetry fix a class of initial 
data for spatial components 

k = -^^-^e.ak'^fir) . (166) 
They contain only one function /(r). The classical equation for this function takes the form 

^f(*.K(*«)=0 ^ d^+ ^^^^^ 
We can see three solutions of this equation. 

/r = 0, /r^ = ±l (r^O). (168) 

^ The author is grateful to V.N. Gribov for the discussion of the problem of instantons during a visit in 
Budapest, May 1996. 



The first solution corresponds to the naive unstable perturbation theory with the assymptotic 
freedom formula. 

Two nontrivial solutions are well-known. They are the Wu-Yang monopoles applied for the 
construction of physical variables in the current literature As it was shown in paper ||2^ the 



Wu-Yang monopole leads to rising potentials of the instantaneous interaction of the quasiparticle 
current. This interaction rearranges the perturbation series, leads to the gluon constituent mass, 
and removes the assymptotic freedom formula ^ as the origin of unstability. 

The Wu-Yang monopole is a solution of classical equations everywhere besides the origin of 
coordinates r = 0. The corresponding magnetic field is 

Bti^k) = ^ . (169) 

Following Wu and Yang pT|, we consider the whole finite space volume, excluding an e-region 
around the singular point. To remove a singularity at the origin of coordinates and regularize 
its energy, the Wu-Yang monopole is considered as the limit of zero size e for the 

Bogomol'nyi-Prasad-Sommerfield (BPS) monopole [^| 



fWY ^ fBPS 

with the finite energy 



1 - 



esinh(r/e) 



(170) 



/Att 1 
£x[B1{^k)? <B^>= — ^ — . (171) 



In this case, the BPS-regularization of the Wu-Yang monopole is the analogue of the infrared 
regularization in QED by the introduction of the photon mass that also violates the initial 
equations of motion. The size of the BPS monopole is chosen so that the parameter e disappears 
in the infinite volume limit 

^ (172) 



as<B'^>V 

and the vacuum energy-density of the monopole solution < 5^ > is removed by a finite counter- 
term in the Lagrangian 

^ . <B'> 
r = r 

2 ■ 

This vacuum magnetic tension is the crucial difference of the topological degeneration of fields 
in Minkowski space from the topological degeneration of classical vacua of the instantons in the 
Euclidean one. 

The problem is to formulate the Dirac quantization of weak perturbations of the non-Abelian 
fields in the presence of the nonperturbative monopole taking into account the topological 
degeneration of all initial data. 

4.4 Dirac method and Gribov copies 

Instead of artificial equations (^), (^) of the gauge-fixing method @ 

F{A^)=0, F{AI) = Mfu^Q ^ Z^^ = J l[DA^detMFS{F{A))e'^ (173) 



we repeat the Dirac constraint-shell formu\at\on resolving the constraint (|140|) 

6W 







D'iA) A'o = D''\{A)doA'r (174) 



with nonzero initial data 

doA^ = =^ A'r{t,x) = <^f\x) . (175) 

The vacuum magnetostatic field ^f^^ has a zero valueof the winding number (|I53|) = 

and satisfies the classical equations everywhere besides of the small region near the origin of 
coordinates of the size of 

e ~ = (176) 

that disappears in the infinite volume limit. 

The second step is the consideration of the perturbation theory ( |163|) where the constraint 
( p74l ) takes the form 

[d^ ($W)]" = 9o[D'^% ($(°))Af . (177) 

Dirac proposed [jl|] that the time component Aq (the quantization of which contradicts to 
quantum principles) can be removed by gauge transformation, so that the constraint ( |177D 
takes the form 

9o[D'^^($(°Vf^] = . (178) 
We defined the constraint- shell gauge 

[D'^^ (<l>(°))Af ^] = (179) 



as the zero initial data of this constraint ( |178|) 



The topological degeneration of initial data means that not only classical vacua but also all 
fields Af^ = + Af^ in the gauge ( p.79| ) are degenerated 

ij") = i;W(f)(if) + di)v^''\xy^ , v^^'^x) = exp[n$o(^)] • (180) 

The winding number functional (|153|) after the transformation ( |162| ) takes the form 

= )] + ^^in) + ^ld'x e^^' TrmAfL-)] , (181) 



where Af{n) = n is given by eq. 

The constraint-shell gauge ( |179|) keeps its form in each topological class 

A"'(4"^)4^"^' = 0, (182) 
if the phase $o(a^) satisfies the equation of the Gribov ambiguity 

[DK^f^)r< = . (183) 

In this case, the topological degeneration means the existence of the manifold of the Gribov 
copies of the constraint-shell gauge ( p.79| ). One can show |^ that the Gribov equation ( |183|) 
together with 

X[$(")]=n (184) 



are compatible with the unique solution of classical equations. It is just the Wu-Yang monopole 
considered before. The nontrivial solution of the equation for the Gribov phase (|183|) in this 
case is well-known 



e 



tanh(r/e) r 



(185) 



it is the Bogomornyi-Prasad-Sommerfield (BPS) monopole 

Thus, instead of the topological degenerated classical vacuum for the instanton calculation 
(that is in the physically unattainable region), we have the topological degenerated Wu-Yang 
monopole ( p.8Cl| ) 

■= v^'^\x)[^f^ + di]v^''\x)-^ , i;(")(f) = exp[n$o(f)] , (186) 
and the topological degenerated multipoles 

A. :=t;(")(f)4. v^''\x) . (187) 

The Gribov copies are evidence of a zero mode in the left hand-side of both the constraints 
O and (III). 

Z}2($(o))]"a^ = . (188) 

A nontrivial solution of this equation 

Al{t,x) = Nmi{x) (189) 

can be removed from the local equations of motion by the gauge transformation (a la Dirac of 
1927) to convert the fields into the Dirac variables 

if) = exp[iV(t)l>o(x)][if) + di] exp[-iV(t)l>o(f)] . (190) 

But this solution (|189| ) cannot be removed from the constraint-shell action W* = J dtN'^I /2+... 
and from the winding number = N + X[A''^^]. Finally we obtained the Feynman path 

integral 

Zf = J DNl[[DEf'''>DAf^]e'^* (191) 

i,c 

that does not coincide with any artificial gauge ( |173|) . 

We consider the derivation of this integral ( |191| ) in details in the following. 



4.5 Topological dynamics 

The repetition of the Dirac definition of the observable variables in QED allowed us to determine 
the vacuum fields and phase of their topological degeneration in the form of the Gribov copies 
of the constraint-shell gauge. 

The degeneration of initial data is the evidence of the zero mode of the Gauss law con- 
straint. In the lowest order of the considered perturbation theory, this constraint ( p.88|) has the 
solution ( |189| ) with the nontrivial vacuum electric field 

F^, = mD'>^i<^Pmx) . (192) 



We call the new variable N(t) the winding number variable as the vacuum Chern-Simons func- 
tional is equal to a difference of the in and out values of this variable 



2 i'out 



u[Ao, $(°)] = ^ J dtj d'xF^^F'^^'^ = ^J d^xF^,B\{^^'^)[N{tout) - iV(t,„)] (193) 
The winding number functional admits its generalization to noninteger degrees of a map p2| 



X[$W] = iV, N^n, (l>f^ = e^*°[<lf +ai]e-^*«) . (194) 

Thus, we can identify the global variable N{t) with the winding number degree of freedom in 
the Minkowskian space described by the action 

W^ = J d'x-{FS,r = J dt^ , (195) 



where the functional 



V 



does not contribute in local equations of motion. The topological degeneration of all fields 
converts into the degeneration of only one global topological variable N{t) with respect to a 
shift of this variable on integers: {N =^ N + n, n = ±1, ±2, ...). Thus, the topological variable 
is a free rotator with the instanton-type wave function (|16C1|) in the Minkowskian space-time 

= exp {zP/viV} , Pn = NI = 2nk + 9 , (197) 

where A; is a number of the Brilloin zone, and 9 is the 6'-angle. In the contrast to the instanton 
wave function (|16C1|) the spectrum of the topological momentum is real and belongs to physical 
values. Finally, equations (|195|) and ( |197|) determine the finite spectrum of the global electric 



tension ( |T^ 

F^, = iV[A($(°Vo]' = (1^ + ^) BK'^^'^) ■ (198) 

It is analogue of the Coleman spectrum of the electric tension in the QED(^i^i^ [^. The appli- 
cation of the Dirac quantization to the 1-dimensional electrodynamics (5-E'-D(i+i) in paper 
demonstrates the universality of the Dirac variables and their adequacy to the description of 
topological dynamics with a nontrivial homotopy group. 

4.6 Zero mode of Gauss Law and Dirac variables 

The constraint-shell theory is obtained by the explicit resolution of the Gauss law constraint 

^ \d'^{A)Y^ A'q = D''\{A)doA'r (199) 



and next in dealing with the initial action on surface of these solutions 



W* = W[A,],^ (200) 



The result of similar solution in QED was the electrostatics and the Coulomb - like atoms. In 
the non-Abelian case, the topological degeneration in the form of the Gribov copies means that 
a general solution of the Gauss law constraint ( |199|) contains the zero mode Z. A general 



solution of the inhomogeneous equation ( |199|) is a sum of the zero-mode solution Z"- of the 
homogeneous equation 

{D'^{A)y^Z^ = , (201) 
and a particular solution Aq of the inhomogeneous one 

A^ = Z'' + A'^ . (202) 

The zero-mode Z^ at the spatial infinity has been represented in the form of sum of the 
product of a new topological variable N(t) and the Gribov phase ^o{x) and weak multipole 
corrections 

Z{t 

) 3^) I asymptotics 

N{t)%{x)+0{-^) , (/> 0) (203) 

In this case, the single one-parametric variable N{t) reproduces the topological degeneration 
of all field variables, if the Dirac variables are defined by the gauge transformations 

= Uz{Z + do)U2' , (204) 

A* = u^iAf^ + d,)U^\ Af)=$f)+Af) 

where the spatial asymptotic of Uz is 

t 

Uz = Texp[ j rft'i(t',f)]|asymptotics = exp [Ar(t) l-^ (f )] . (205) 

The topological degeneration of all fields converts into the degeneration of only one global 
topological variable A^(t) with respect to a shift of this variable on integers: (A^ =^ N + n,n = 
±1,±2,...). 

4.7 Constraining with the zero mode 

Let us formulate an equivalent unconstrained system for the YM theory in the monopole class 
of functions in the presence the zero mode Z^ of the Gauss law constraint 

A-, = Z'^ + = -Dt\A)Z' + F,\ ( {D\A)rz'' = ) . (206) 

To obtain the constraint-shell action 

W^YAf (constraint) = WymI^] + Wym[F] , (207) 
we use the evident decomposition 

= {-DZ + Ff = {DZf - 2FDZ + {Ff = d{Z{DZ)) - 2d{ZF) + {Ff (208) 
and the Gauss Eqs. DF = and D^Z = which show that the zero mode part Wym of the 



constraint-shell action (|207| ) is the sum of two surface integrals 



Wym[2] = j dtj d'x[^d,iZ''DfiA)Z') - 9,(Fo",Z'^)] = W° + W , (209) 
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where the first one W° is the kinetic term and the second one W describes the coupling of the 
zero-mode to the local excitations. These surface terms are determined by the asymptotic of the 
fields iZ'',A^) at spatial infinity ( PU3|) , (|TBg) which we denoted by (iV(t)<l>g(f), <l>^(f)). The 



fluctuations Fq\ belongs to the class of multipoles. Since the surface integral over monopole- 
multipole couplings vanishes, the fluctuation part of the second term obviously drops out. The 
substitution of the solution with the asymptotic ( |203| ) into the first surface term Eq. ( [209| ) leads 
to the zero-mode action ( |195|) . 

The action for the equivalent unconstrained system of the local excitations, 

WymIF] = J d'x [e^, . Af^ - \ [El + + [Df , (210) 

is obtained in terms of variables with zero degree of map 



-~(o) 

TT T:^\^ITT~L a TT I A\^) I Q \TT~l- A \^ ) I J. Zf\ I / i 

Ok 

AO) 



UzF^l^U^' , A = UziAf'^ + d,)U^' , Af'\t,x)=¥f\x)+A {t,x) . (211) 



by decomposing the electrical components of the field strength tensor Fgj into transverse E°- 
and longitudinal F^-^ = -Df ($W)cr^ parts, so that 

°) = - (<l>(°))a''; ($(°))E^ = . (212) 

Here the function is determined from the Gauss equation 

((Z}2($(o)))a(> + (^e^'^-if ($(0))) = -ge'^^'Af^E^ . (213) 

Due to gauge-invariance the dependence of the action for local ecxitations on the zero mode 
disappears, and we got the ordinary generalization of the Coulomb gauge 0, ||] in the presence 
of the Wu-Yang monopole. 

4.8 Feynman path integral 

The Feynman path integral over the independent variables includes the integration over the 
topological variable N{t) 

Zf[J]= fUdNit)Z[J^] , (214) 

-I- 



where 

ZiJ""] = /n {n l^^^^^^^^l exp^{WyMiZ) + WyM{A^^^) + S[J^]} . (215) 



t,x \a=l 



As we have seen above, the functionals W, S are given in terms of the variables which contain 
the nonperturbative phase factors U = Uz ( |205| ) of the topological degeneration of initial data. 
These factors disappear in the action W, but not in the source 

S[J^] = J d^xJ^Al A, = U{Af)U-^ (216) 

what reflects the fact of the topological degeneration of the physical fields. 

The constraint-shell formulation distinguishes a bare "gluon" as a weak deviation of the 
monopole with the index [n = 0), and an observable (physical) "gluon" averaged over the 
topological degeneration (i.e., Gribov's copies) [|l7j 

A^'^y^ = lim ^ ^f: ^^"^ (^) - Ko ; (217) 

n= — L 

whereas, in QED the constraint-shell field is a transversal photon. 



4.9 Rising potential induced by monopole 

We can calculate the instantaneous the Green function 

{D\^^'^^)y\x)G''%x, y) = d^'d^x - y) . (218) 
In the presence of the Wu-Yang monopole we have 

and na{x) = Xa/r; r = \x\ . Let us decompose G"^ into a complete set of orthogonal vectors 
in color space 

C^^ix^y) = [n'^ix)n\y)Vo{z) + e^(x)e^(y)\/i(z)]; {z = \x - y\) . 

a=l,2 

Substituting the latter into the first equation, we get 

f^K + -:fK-4^n = o n = 0, 1. 

az"^ z az z'^ 
The general solution for the last equation is 

=rf„|f-y|'" +c„|x-y|'" , n = 0, 1, (219) 

where d^, Cn are constants, and can be found as roots of the equation (/")^ + /" = n , 

It is easy to see that for n = we get the Coulomb-type potential = —l/An, 

;? = -i^ = -i; i5 = ^ii^ = o, (221) 

and for = 1 the 'golden section' potential with 

= _i±}/? ^ _1.618 ; 4 = Zi±^ ^ 0.618 . (222) 

The last potential (in the contrast with the Coulomb-type one) can leads to rearrangement of 
the naive perturbation series of the type of the spontaneous chiral symmetry breaking. This 
potential can be considered as the origin of the 'hadronization' of quarks and gluons in QCD 



4.10 FP path integral 

Thus, we can say that the Dirac variables with the topological degeneration of initial states in a 
non-Abelian theory determine the physical origin of hadronization and confinement as nonlocal 
monopole effects. The Dirac variables distinguish a unque gauge. In QED, it is the Coulomb 
gauge; whereas, in YM theory, it is covariant generalization the Coulomb gauge in the presence 
of a monopole. 
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If we pass to another gauges of physical sources on the level of the FP integral in relativistic 
gauges, all the monopole effects of the degeneration and rising potential can be lost (as the 
Coulomb potential is lost in QED in relativistic invariant gauges). Recall that to prove the 
equivalence of the Feynman integral to the Faddeev-Popov integral in an arbitrary gauge, we 
change variables and concentrate all monopole effects in the phase factors before the physical 
sources. The change of the sources removes all these effects. 

The change of the sources was possible in the Abelian theory only for the scattering ampli- 
tudes when all particle-like excitations of the fields are on their mass-shell. However, for the 
cases of nonlocal bound states and other phenomena where these fields are off their mass-shell 
the Faddeev theorem of equivalence of different "gauges" is not valid. 

For the non-Abelian theory all states become nonlocal, and the range of validity of the 
equivalence theorem is equal to zero. 

5 Dirac variables in QCD: hadronization and confine- 
ment 

5.1 Action and first principles of QCD 

In the previous Sections we have shown that the Dirac variables of gauge theories is the best 
way to describe both the topological degeneration of initial data and the relativistic covariant 
bound states. 

Here we would like to demonstrate the adequacy of the Dirac variables to the physical 
phenomena of hadronization and confinement in QCD. 

QCD was proposed [|80|] as the non-Abelian SUc{S) theory with the action functional 

W = J (i'x |^(GS,' - Bf) + ij[zr{d, + 4) - m]^^J^ , (223) 

where ip and ip are the fermionic quark fields, 

GS, = doA^ - Df{A)Al , = e.,fc (d,At + ^-r"'A]Ai) (224) 

are non-Abelian electric and magnetic fields, and D'^^^A) := + gf^'^^Al is the covariant 
derivative. 

The action (^2^) is invariant with respect to relativistic and gauge transformations u{t,x) 
A1 ■=u{t,x)[A + di)u-\t,x), i)"" ■.= u{t,x)il) , (225) 
where A^ = g^A^ . 

5.2 Topological degeneration of QCD 

The group of stationary gauge transformations in QCD represents the group of three-dimensional 
paths u{t,x) = v{x) lying on the three-dimensional space of the S'f/c(3)-manifold with the 
homotopy group 7r(3)(5'f/c(3)) = Z. Thus, a difference between the Abelian and non-Abelian 
theories is the difference between the homotopy groups of maps of three-dimensional coordinate 
space onto the gauge manifolds f/(l), and SU{3)c respectively. 
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The whole group of stationary gauge transformations is split into topological classes marked 
by the integer number n (the degree of the map) which counts how many times a three- 
dimensional path turns around the S'f/(3)-manifold when the coordinate Xi runs over the space 
where it is defined. The stationary transformations t>"(x) with n = are called the small ones; 
and those with n ^ 



(226) 



the large ones. 

The degree of a map 



U[n] 



247r2 



n 



(227) 



as the condition for normalization means that the large transformations are given in the class of 
functions with the spatial asymptotic 0{l/r). Such a function L" ( |226| ) is given by 



)(f) = exp(nl>o(f)), $0 



where the antisymmetric SU(3) matrices are denoted by 



-^vr^/o(r) 



(228) 



:— A^, \\ :— A^, A^ :— A''' , 
and r = \x\. The function /o(?") satisfies the boundary conditions (|157|) 

/o(0) = 0, /o(oo) = 1 . 



(229) 



The functions L" belong to monopole-type class of functions. The topological degeneration of 
QCD means that initial data of the transformed physical fields Ai in ( |226| ) should be also given 
in the same monopole class of functions. 

We have seen in the previous Section that, in the monopole class of functions in the sector 
with a zero value of the winding number functional, there is the gauge Higgs effect. This effect 
leads to a stable physical vacuum in the form of the 5't/(3) BPS monopole 



($ 



(O)nBPS 



X 



BPS 



2 

where e is defined by the monopole energy 

in 1 



ft 



BPS 



1 - 



esinh(r/e) 



(230) 



d'x[Bn^k)y 



y <B^> 



a. 



a 



47r' 



V 



(TX 



(231) 



The average vacuum magnetic tension $ 



< B' >-- 



(0), 



O(-) with a zero winding number and 



Jd'^x te($(°)) 



/ d^x 



7^0 



(232) 



is the main difference of the consiramt-s/ie// formulation of non-Abelian theory in the Minkowski 
space-time from the instanton calculation with functions interpolating between two classical 
vacua with zero energy that belong to nonphysical values of spectrum in the quantum theory. 
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Any field is a sum 

A^\t,x) = <^f\x) + Af\t,x) , (233) 
where Af'^ is a weak perturbative part with the asymptotic at the spatial infinity 

A^°^(t,f)|asymptotics = O(^) (/ > l) . (234) 

The Dirac variables distinguish the constraint-shell gauge 

Dl^i^f^)Af\to,x)= 0. (235) 

The topological degeneration 

= exp[n%{x)] [Af^ + dk) exp[-n%{x)] (236) 
keeps the constraint-shell gauge ( |235| ), if there is a nontrivial solution of the Gribov equation 

A'($f)<fo = 0. (237) 
The nontrivial solution of the Gribov equation exists in the form of the SU (3) BPS monopole 

1 e 



tanh(r/e) r 



(238) 



This solution has the boundary condition ( |229|) of a phase of the topological transforma- 
tions (|2|). 

We get the vacuum as an exact solution in a form of the Wu-Yang monopole in the 
limit of an infinite volume e — > at the end of the calculation of spectra and matrix elements 

lim = 1 . (239) 

5.3 Superfluid dynamics of gluonic liquid 

We have seen above (in Section 4) the constraint-shell action can contains the superfluid quan- 
tum motion of the gauge phases considered as a zero mode {Z) of the Gauss law constraint 

^ = ^ {D\A)rAo' = DriA)doA'r + , (240) 
dAQ 

where = gip^'jf^ip is the quark current. In lowest order of the perturbation theory this 
constraint takes the form 

(D|($))»'^Ao'^ = Dr($)ao4^ , (241) 



In the considered case of BPS-monopole ( |230|) , the Gribov copies mean that there is the zero 



mode of the covariant Laplace operator in the monopole field 

{DY\^k^'^)i^o^'^)\x) = . (242) 



This is a perturbative form of a zero mode of the Gauss law constraint ( |24C1| ) [|I7| , as the 
solution of the homogeneous equation ( |24C1| ) 



{D\A)rZ'' = , (243) 



with the asymptotics at the space infinity 



i(t,x)|asymptotics = iV(t)<lo(x) , (244) 

where N(t) is the global variable of an excitation of the gluon system as a whole. 

From the mathematical point of view, the zero mode means that the general solution of the 
inhomogeneous equation ( [24C1| ) for the time-like component Aq is a sum of the homogeneous 
equation ( |243| ) and a particular solution Aq of the inhomogeneous one ( |240| ) : Aq = + Aq . 

The zero mode of the Gauss constraint and the topological variable N(t) allow us to describe 
the topological degeneration of all fields by the non-Abelian generalization of the Dirac dressed 
variables in QED 

= Uz{Z + doWz' , A* = Uz{A^'^ + di)Uz\ r = Uzi^^'^ , (245) 

where the spatial asymptotic of Uz is 

t 

Uz = Texp[ J dt'Z{t\ f)]|asymptotic. = exp [iV(t) $0 (x)] = f/(f ) , (246) 

and A^^^ = + A^^\-ip^^^ are the degeneration free variables. In this case, the topological 
degeneration of all color fields converts into the degeneration of only one global topological 
variable iV(t) with respect to a shift of this variable on integers: (A^ =^ N+n, n = ±1, ±2, ...). 

One can check that the Pontryagin index for the Dirac variables ( |245| ) with the asymp- 
totic (|234|) , ( [244| ) , ( |246| ) is determined by only the difference of the final and initial values of 
the topological variable 

2 tout 

i^[A*] = ^ I dtj d'xGl^^G'^'^'' = N{t,^,) - N{U^) . (247) 

tin 

Thus, we can identify the global variable N(t) with the winding number degree of freedom in 
the Minkowski space. 

The dynamics of physical variables including the topological one is determined by the constra- 
int-shell action of an equivalent unconstrained system as a sum of the zero mode part, and the 
monopole and perturbative ones 

W*,,, = Wcauss-shell = Wz[N] + W^on[^,] + Wioc[A] . (248) 

The action for an equivalent unconstrained system ( |248|) in the gauge 



Dl^{^)Al = . (249) 



with a monopole and a zero mode has been obtained in the paper [^. This action contains 
the dynamics of the topological variable in the form of a free rotator 

. ,.iV2j . /■ 47r 
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//V / r 4-Tr 

dt^; I = J d^Dn^km' = -(2vr)^e , (250) 

V ^ 



where e is a size of the BPS monopole considered as a parameter of the infrared regularization 
which disappears in the infinite volume limit. The dependence of e on volume can be chosen as 
e ~ V^^, so that the density of energy was finite. 



The perturbation theory in the sector of local excitations PVioc[^] is based on the Green 
function in the form of the inverse differential operator of the Gauss law 



[D\<!>)]'''V'\x, y) = -6\x - y)6'"' (251) 

which is the non-Abelian generalization of the Coulomb potential. As it has been shown in the 
previous Section the non-Abelian Green function in the field of the Wu-Yang monopole is the sum 
of a Coulomb-type potential and a rising one [|^. This means that the instantaneous quark- 



quark interaction leads to spontaneous chiral symmetry breaking Goldstone mesonic 



bound states p5[|, rearrangement of perturbation series for gluons, glueballs and the Gribov 



modification of the asymptotic freedom formula. If we choose a time-axis l^^^ along the total 
momentum of bound states (this choice is compatible with the experience of QED in 
the description of instantaneous bound states), we get the bilocal generalization of the chiral 



Lagrangian-type mesonic interactions [^. In this case, the U(l) anomalous interaction of 770- 
meson with the topological variable lead to additional mass of this isoscalar meson ||2^. All 
these effects can be described by the Feynman integral. 

5.4 Feynman path integral 

The Feynman path integral for the obtained unconstrained system in the class of functions of 
the topological transformations takes the form (see p2| ) 

, c=8 



ZF[l^°\r*] = [ DN{t) [ l[[d^A'*d^E''*]l[[dij*d^*] 

•' •' c=l 

xexp{zW*^o)[N,A*,E*,^p*,ij*]+zS*} , (252) 

where W* is given by Eq. (|248|) , 

S* = J dS[r; ■ A"* + s*ij* + ij*s*] , (253) 

and J'^*,s*,s* are physical sources. 

The action of the local variables does not depend on the Gribov phases 

WUA*) + ... = WUA^''^) + ... . (254) 

Whereas the physical sources of the Dirac dressed variables A*,^* are topologically degener- 
ated. The nonperturbative phase factors of the topological degeneration can lead to a complete 
destructive interference of color amplitudes [|17], ^ due to averaging over all parameters of the 
degenerations. 

The perturbation theory is compatible with only one gauge. It is the Coulomb-type gauge 
in the monopole field [Q. 

The change of variables A* concentrates all monopole physics in the form of the non-Abelian 
Dirac factor 

U{A) = f/^ |l + -^^D,{^)A, + ...| (255) 

before the " physical" sources. The change of " physical" sources A* J* =^ (that is called the 
transition to another gauge) removes all monopole physics, including confinement and hadroniza- 
tion, like the similar change in QED removes all electrostatic phenomena from the FP integral 
in the relativistic gauges. 



5.5 A free rotator: topological confinement 



We have seen in Section 4, that the topology can be the origin of color confinement as complete 
destructive interference of the phase factors of the topological degeneration of initial data. 

The mechanical analogue of the topological degeneration of initial data is the free rotator 
N(t) with the action of free particle 

W{NouuN,n\ti) = ] dt^I, p = NI, H^ = Yi (256) 



given on a ring where the points N{t) + n [n\s integer) are physically equivalent. Instead of a 
initial date N{t = 0) = in the mechanics in the space with the trivial topology, the observer 
of the rotator has the manifold of initial data A^(")(t = 0) = A^i„ + n; n = 0, ±1, ±2, .... 

An observer does not know where is the rotator. It can be at points Nin, Nin ± 1, A^m ± 
2, Nin ± 3, .... Therefore, he should to average a wave function 

^(N) = e'P^ . 

over all values of the topological degeneration with the 6'-angle measure exp(i^n). In the result 
we obtain the wave function 

1 n=+L 

^(iV)observabie = Hm — e''"^(A^ + ^) = exp{t{27Tk + e)N} , (257) 

L^oo ZL — T 
n=—L 

where k is integer. In the opposite case p ^ 2Txk + 6 the corresponding wave function (i.e., the 
probability amplitude) disappears \&(iV)observabie = due to the complete destructive interfer- 
ence. 

The consequence of this topological degeneration is that a part of values of the momentum 
spectrum becomes unobservable in the comparison with the trivial topology. 

This fact can be treated as confinement of those values which do not coincide with the 
discreet ones 

Pk = 2-nk + e, Q<e <Ti . (258) 

The observable spectrum follows also from the constraint of the equivalence of the point N and 

N + 1 

m{N) = e'^-^{N + 1), ^(A^) = e'P^ . (259) 

In the result we obtain the spectral decomposition of the Green function of the free rotator ( |256|) 
(as the probability amplitude of transition from the point Nin to Nout) over the observable values 
of spectrum ( |258| ) 

fc = + oo 

G{Nout, Nin\ti) =< Nout\exp{-iHti)\Nin >= 7^ exp 

k=—oo 

Using the connection with the Jacobian theta-functions [^] 

k=+oo 

e3(Z|r) = V exp[i7ik^T + 2ikZ] = {-ir)-^^^ exp[—]e3 

fe=— 00 



-i^h + ipk{Nout - Nin) ■ 

(260) 
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we can represent expression ( [26C1| ) in a form of the sum over all paths 



/ J n=+oo 

G{Nout, Ninlh) = J —— exp[i9n]exp[+iW{Nout,Nin + n\ti)] , (261) 

where 

UK AT J. AT (Nout + n-Ninfl 



2ti 



is the rotator action ( |256| ). 



5.6 Confinement as destructive interference 

The similar topological confinement as a complete destructive interference of phase factor of the 
topological degeneration (i.e., a pure quantum effect) can be in the "classical non-Abelian field 
theory". Recall that at the time of the first paper of Dirac ^ (1927) the so called "classical 
relativistic field theory" were revealed in the papers of Schrodinger, Fock, Klein, Weyl [^, ^ 



as a type of relativistic quantum mechanics, i.e., the result of the primary quantization. The 
phase of the gauge transformations was introduced by Weyl as a pure quantum quantity. 

The free rotator shows us that the topological degeneracy can be removed, if all Green 
functions are averaged overvalues of the topological variable and all possible angles of orientation 
of the monopole unit vector {n = x/r) ( |164| ) in the group space (instead of the instanton 
averaging over interpolations between different vacua). 

Averaging over all parameters of the degenerations can lead to a complete destructive inter- 
ference of all color amplitudes |]T^, [7T|. In this case, only colorless (" hadron" ) states have to 
form a complete set of physical states. Using the example of a free rotator, we have seen that 
the disappearance of a part of physical states due to the topological degeneration (confinement) 
does not violate the composition law for Green functions 

Gij{ti,t3) = 'YGih(ti,t2)Ghj{t2,t-i) (262) 

h 

defined as the amplitude of the probability to find a system with the Hamiltonian H '\n a state j 
at the time t^, if at the time ti this system was in a state i, where {i,j) belongs to a complete 
set of all states {h}: 

Gij{ti,t^) =< i\exp{-i J H)\j > 

tl 

The particular case of this composition law ( p62| ) is the unitarity of S-matrix 

SS+ = 1 ^ ^ < i\S\h >< h\S-^\j >=< i\j > 

h 

known as the law of probability conservation for S-matrix elements {S = I + iT) 

J2 < i\T\h >< h\T*\j >= 2lm < i\T\j > . (263) 

h 

The left side of this law is the analogue of the spectral series of the free rotator (|260|) . The 
destructive interference keeps only colorless "hadron" states. Whereas, the right side of this law 
far from resonances can presented in a form of the perturbation series over the Feynman diagrams 



47 



that follow from the Hamiltonian. Due to gauge invariance H[A^"'\ q^'"'^] = H[A^^\ q^^^] this 
Hamiltonian does not depend on the Gribov phase factors and it contains the perturbation series 
in terms of only the zero - map fields (i.e., in terms of constituent color particles) that can be 
identified with the Feynman partons. The Feynman path integral as the generation functional 
of this perturbation series is the analogue of the sum over all path of the free rotator ( |261j ). 

Therefore, confinement in the spirit of the complete destructive interference of color ampli- 
tudes |]T^, |T7|, |7I[| and the law of probability conservation for S-matrix elements ( |263| ) leads to 



Feynman quark-hadron duality, that is foundation of all the parton model and QCD appli- 
cation [^. The quark-parton duality gives the method of direct experimental measurement of 
the quark and gluon quantum numbers from the deep inelastic scattering cross-section [^. For 
example, according to Particle Data Group the ratio of the sum of the probabilities of r-decay 
hadron modes to the probability of r-decay muon mode is 



h 



3.3 ±0.3 . 



Wr- 



This is the left-hand side of Eq. ( |263|) normalized to the value of the lepton mode probability 
of r-decay. On the right-hand side of Eq. (|263|) we have the ratio of the imaginary part of the 
sum of quark-gluon diagrams (in terms of constituent fields free from the Gribov phase factors) 
to the one of the lepton diagram. In the lowest order of QCD perturbation on the right-hand 
side we get the number of colors and, therefore, 

3.3 ± 0.3 = . 

Thus in the constraint-shell QCD we can understand not only "why we do not see quarks", 
but also "why we can measure their quantum numbers". This mechanism of confinement due 
to the quantum interference of phase factors (revealed by the explicit resolving the Gauss law 
constraint [l^, |7l|, |l9l) disappears after the change of "physical" sources A* J* =^ AJ that is 
called the transition to another gauge in the gauge-fixing method. 



5.7 Relativistic equations for gluonic systems 
The Wu-Yang monopole 



does not change qualitatively the spectrum of gluon excitations in the comparison with the naive 
perturbation theory 

= dlA^ ■ 



There are only the spin-color mixing [^. Here we omit this mixing and consider mainly the 
instantaneous current-current interaction described by the QCD Hamiltonian ^ 



H- 



= J d^lmr + m'] + 

Here V^(x) denotes the sum of rising potentials induced by the Wu-Yang monopole. 
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Let us represent the gluon fields as Bogolubov expansion in creation and annihilation oper- 
ators 



^r(x) = [ I \ faW^(p)er + a(-)''(p)er]e-P-, (264) 

where the function 0(p) is to be calculated from the diagonalization condition of the Hamilto- 
nian ( [264| ) with respect to the operators a'^'^\a^^\ 
In the normal form, the gluon Hamiltonian reads as 

cik[ ( aW(-k)aW(k)+aL7)(k)aL;)(-k)) ^ + 

+ aWai;)c.--(0) + O(a^)]. (265) 

Here C(0) and ct;(</)) are some defined functions the dependence on of which is determined 
from the diagonalization of the Hamiltonian ( [265| ). 

The diagonalization condition for the Hamiltonian ( |265|) means that the coefficients C(0) 
vanish 

C(0) = (266) 

For the solutions ( |266|) the function 07(0) defines the gluon energy spectrum in the same way 
as the Schwinger-Dyson equation defines the energy spectrum for the quarks. 

Let us illustrate this scheme by calculating the one-particle energy of the gluon and its 
bound states for the simplest example of the theory ( [2b4| ) where the operator for the potential 
is substituted by an effective potential. This means we consider the sum of the free Hamiltonian 

and the Hamiltonian of potential interaction between colour gluon currents 



(pi)0(qi) 



(27r)3 \\ 0(p2)0(q2) 
• _ a(-)(-l)][aW(2) + a^-\-2)] ■ 

■ [aW(-l') - + a^"^(-2')]- 

Here, the short - hand notation (±pi, 61, ii) = (±1), (±qi, ci, ji) = (±2), 

(±P2,62,^2) = (il'); (='=q25 C2, j2) = (=1=2') has been used. In our case, the coefficients C and 

LU are given by the expressions 



C"^"^(0) = [ ^'^'' ^ iM\W _ ^;,^.]e^i(k)ef (-k) = 0, 



o^"^"^(0) = [^^^^^^^^^^^ + ^^^]er(k)ef(k) (267) 



4Q 



with 

= f^/ (^Hk - q)0(q)(iy - ^ii?,), 

Since two gluons form the simplest bound state, by analogue with the mesons, we introduce 
a glueball creation operator 

= E /^k [ X(+t)(k)e7Hk)ef (-k)aW''(k)aW\-k) + 
b ■' 

+ ^fc;nk)er(k)e7(-k)a(7)^(k)a(;)^(-k)] 

and a "coherent" vacuum 

|0»,-exp{E / rfkirfk2a(ki,k2)(a(+)^(ki)aW^(k2))- 

c 

• «^''(k0aW''(k2))}|0 > . 
Then, the Schrodinger equation for eigenvalues of the Hamiltonian operator 

a « 0|G'ii'G'+|0 Mg a « OIGG+IO »«, (269) 

is equivalent to the equation for the glueball wave function 

(2^(k) - Mg)x\;'-\\.) 



(2u;(k) + MG)x\--\]<i) 

^'^:4,{(W^+(q|k))2x(r-)(q) - {W-{c^\\.)fxl,~\c^)} 



with 



iy±(q|k) 



(270) 



(271) 



W(q) = /^Hk-q)/(q). (272) 

Furthermore, one can write also the equation for the wave function of a antisymmetric three - 
gluon state '0+++('0+ , -0-++, V'-+-) with eigenvalues Msg 

[uj{l) + a;(2)+a;(3)-MsG]^+++(l,2,3) = 

N ^ 

= ^{/i,2(a;MV'+++(i,2,3)+(^ria;22V^__+(l,2,3))+ ... }. (273) 



For an estimate of the solution to equations ( [27C1| ), ( [272| ) we make use of the separable 
approximation for the potential 

with the parameters L = l.QGeV, figcD = O-SbGeV. Equation ( P69|) take the form 



1 + ^/^*^^. (275) 



m2 



^(p) = V^A/p^-m2 
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and have the solutions ~ O.SGeV, ^ 1.18. For the scalar glueball mass as an eigenvalue 
of equations ( |270|) , ( |272| ) in the separable approximation 



[2Vz^k^ + ml - Mg]X 

Nc dq 



I -7^sli^-Y)W{k\q) + 2{X + Y)], 



[2v^.^k2 + ml + Mg]X = 

' J 7A[(>^-^)W^(k|q) + 2(X + r)], (276) 



W{k\q) : 



q2 + 



mi 



one obtains the value Mq — l.GGeV. 

The Green function for the transversal gluon corresponding to equations ( |265|) and ( |266|) is 
given by 

Ai(go, q) = -(^ij - Qi^Qj)- (2T7) 

0(q) go (q) q 

From its meaning the quantity ^ = Z(q) can be called the infrared renormalization constant of 
the wave function. 

The Green functions for the constituent quarks and gluons are elements of a new quasi particle 
perturbation theory in terms of which all matrix elements are calculated including the " running" 
coupling constant |p9 |. 



The phenomenon of dimensional transmutation appearing in the " running" coupling constant 
should be investigated in accordance with the logic of quantum theory at the stage of defining the 
quark and gluon energy spectrum and their one-particle Green functions. We have seen before 
that in QCD like in QED the infrared parameters appear as nontrivial boundary conditions of 
the equations of motion. 

The appearance of constituent masses for quarks and gluons does of course influence the 
determination of the "running" coupling constant which in the new theory cannot have any 
singularities in the whole Euclidean region of the transferred momenta, among them also at 
q^ = |3l|, |50l]. Thus, the rising potential rearranges the perturbation theory and makes it 
stable. 



5.8 Hadronization and chiral Lagrangian limit 

One of our supposition was that the main role in the formation of instantaneous bound state in 
gauge theories plays the "instantaneous" interaction. In QCD this interaction 

Wins. = -I J d'xd'yijfMi^^^\{x)JC2S^^^^ (278) 

is determined by the Green function ( |251D as the inverse differential operator of the Gauss law 
constraint G"^ in the field of the Wu-Yang monopole 



which is the non-Abelian generalization of the Coulomb potential. We have seen above in 
previous Section that the non-Abelian Green function in the field of the Wu-Yang monopole is 
the sum of a Coulomb-type potential and a rising one 



The procedure of the hadronization, as we have seen in Section 3 consists in the application 
of the Legendre transformation with bilocal field Ai{x,y) = Ai{z\X) 

^ / d^xd'^y{tlj{y)^{x)))C{x,y){tlj{x)^{y)) = 

= '1 I dVyMix,y)IC~\x,y)Mix,y)+ (280) 
+ / d'^xd^y{ij{x)ij{y))M{x,y) 

where JC^^ is the inverse of the kernel ( |278| ). In the short- hand notation 

d'^xilj{x){ip{A^) - m^)tfj{x) = f d^xd^ytlj{y)ij{x){ip{A^) - m^)6{x - y) 



= {^|J^P,-G-\A^)) , 
d'^xd'^y{ijj{x)ip{y))M{x,y) = {ijjijj^M) 

the QCD action is represented in the form 

WQCD = Wg{A^) + W[^,M\A^] , (281) 

where 

W[^p, M \A^] = (#, i-G-\A^) +M))-^iM, IC-^M) , (282) 

and Wg{A-^) is the action for transverse gluons. After quantization ( or integration ) over quark 
fields and normal ordering, the action takes the form 

W[M\A^] = -^{M,IC-'M)+iJrln [Gm{A^)] (283) 

where the Green function Gj^iA-^) is defined by equation 

[G{A^)-^ - M]Gm{A^) = I (284) 



or 



d^y[{ip{A^) - m'^)6\x - y) - M{x,y)\GM{y, z\A^) = 6\x - z) (285) 

As a result of such quantization, only the the contributions with inner fermionic lines (but 
no scattering and dissociation channel contribution) are included in the effective action since we 
are interested only in the bound states. 

Finally, we got the following result 

Z[. ..]=<! ciA<e^^^eff(^l^")+-> (286) 

here < F(A-^) >^x is the Feynman integral ( |252|) depended on the time- axis r]^ according to 
Markov- Yukawa prescription. 

The determination of a minimum of the action ( |283| ) leads to the SD equation 

^^^^^ = ^ S(A^) = ^}CG^{A^) . (287) 

We denoted the corresponding classical solution for the bilocal field by S. 

The next step is the expansion of the action (|283| ) around the point of minimum Ai = 

? 

Weff{J: + M') = Weff{J:) + [-^M'IC-'M'+'-{G^MT] + 

+ ^E-(G'e-M')", , (288) 

where 

GM{x,y) = j (fzG{x,z)M{z,y) = <^{x,y), 

$2 = tr J d'^xd'^y<^{x,y)<^{y,x), (289) 
$^ = tr J d^xd'^yd*z^{x,y)(^{y, z)^{z,x) . . 
The small fluctuations Ai' we represent as a sum over the complete set of classical solutions F, 

6'Wq{j: + m') 

\m'=o ■ i - . (290) 

The bound state function T{x,y\A-^) satisfies the Bethe- Salpeter equation in the form 

T{x,y\A^) = i)C{x,y) J d^z,d'z2Gj,{x, z,\A^)T{z,Z2\A^)Gj:{z2,y\A^) . (291) 
By the same way we can construct any conventional vertex function 

Tcix,y\A^) = t}Cix,y) [ d^z,d%G^ix, z,\A^)^c{ziZ2\A^)Gj:{z2,y\A^) . (292) 



The 7c is one of the complete set of the Dirac matrices. 

The obtained expression are the exact equations for the description heavy and light meson 
in QCD in which we may to try to separate the radiation corrections (retardations) from bound 



state effects. All radiation corrections are separated on two parts: the "retardation" and the 
"instantaneous" (forming the instantaneous potential). 

We can consider the obtained expressions in the lowest order in the retardation for the 
colourless hadron 

Z = j dMe^^^fi^^^^ , (293) 



and 



1 °° 1 

W[M\ = --N,{M, IC-'M) + iN, , (294) 

^ n=l ^ 

where /C and $ are the colourless operators 

= ^^fiV^^f{z^)6{z . r^)f, . (295) 

In this approach, we get the standard Schwinger- Dyson (SD) and Bethe- Salpeter (BS) 
equations: 

- y) = m°6^^\x - y) + i!C{x,y)Gj:{x - y), (296) 

r = z/C(a;, y) J <fzidh2Gj,{x - zi)T{z^, ^2)^2(^2 - y) (297) 

They describe the spectrum of Dirac particles in bound states and the spectrum of the bound 
states themselves, respectively, like in Sections 3.3.-3-5.. 

The relativistic bound state equations derived from the constraint-shell gauge theory with 
a rising potential contains the pure quantum effect of rearrangement of the ground state with 
spontaneous breakdown of chiral symmetry. This effect leads to the low - energy theory of light 
mesons, in which the pion is considered in two different ways, as a quark - antiquark bound 
state and as a Goldstone particle. 

There is a number of paper (cf. ^ 0, |5^, |58|, ^ and references therein) 



where SD and BS equations are used for the calculation of the mass spectrum of light mesons, 
the constituent quark masses and the meson decay constants. In the papers the potentials are 
determined from the spectroscopy of heavy quarkonia as sums of rising and Coulomb potentials, 
for instance 

Vir) = ^ - ^r^, Vn'^ - 250MeV, as ~ 0.3 . (298) 

r 

Thereby, the heavy quarkonia (m° >> 250MeV) themselves are described by Schrodinger 
equation which, as has been shown, can be derived from BS equation in the limit of large 
masses. In this limit, the effect of spontaneous breakdown of chiral symmetry also disappears, 
and the constituent quark masses are identical to the current ones [^, 



The advantage of such a potential approach compared with all other ones consists in the 
first constructive connection between the fundamental parameters for physics at short distances 
(the parameters of rising and Coulomb potentials and the current quark masses) with those of 
hadron physics for long distances (the pion mass and its weak decay constant F^). 

The shortcomings of this approach were the following: the nonrelativistic formulation (in the 
rest frame) of the bound state, the absence of an relativistic meson interaction Hamiltonian, and 



the open problem of the status of radiative QCD corrections. The first two disadvantages are 
absent in the new relativistic potential model ^ considered here on the basis of the Dirac 
variables and the Markov- Yukawa prescription. This model represents a logical interpretation of 
relativistic atomic physics, i.e. an interpretation of the " atomization" of QED. 

From this point of view the " hadronization" of QCD qualitatively differs only by the short 
- range property of the quark - antiquark interaction potential for light quarkonia. Furthermore, 
the effective action for light mesons must be an action for a chiral Lagrangian The proof 



of the fact that the effective bilocal action leads to a chiral Lagrangian has been performed 
in with the help of the separable approximation which can be used just for short - range 
potentials. For low orbital momenta such potentials can be represented with good accuracy as 
a product of two factors 

<l = 0\V{p~q)^\l = 0>=f{p^)f{q^). 
The underlying bilocal model becomes equivalent to one of versions of the Nambu - Jona - Lasinio 



model ^] with explicit indication of the form factor f (p) for the ultraviolet regularization. 



It is well known |^, that this model leads to chiral Lagrangians. 

The validity of the separable approximation for short - range potentials explains the fact 
of the weak dependence of low - energy physics for light mesons on the form of the potential. 
Therefore, there exists a number of models yielding a satisfactory description of the experimental 
data. 

Here, one should mention also papers dealing with the derivation of nonlinear chiral La- 
grangians from QCD (cf.ref. in ||3^). The essence of those proofs consists in a formal derivation 
of this Lagrangian by means of chiral transformation which are parameterized by the meson field. 
Thereby, in many cases no derivation of the equation for the meson spectrum is given to say 
nothing of its solution. The main aim of these papers is to find the coefficients in higher order 
terms in the expansion of chiral Lagrangians in meson momenta and to establish the description 
of baryons in the form of "skyrmions" [^. All these papers concerning the justification of chiral 
Lagrangian from QCD are not devoted to the determination of essential parameters of the low 
- energy physics ( Fj,,FK,m^, ...) from QCD. 

The relativistic bilocal model for atoms and hadrons as compared with the above - men- 
tioned popular non - relativistic [^, ^ and nonlinear approaches unifies aspects of both 



approaches and gives a constructive generalization of chiral Lagrangians to heavy quarkonium 
physics, i.e. it allows to describe decays of heavy quarkonia into light ones in the framework 
of the bilocal action with a minimal number of parameters, defined in the short - range region 
where the perturbation theory begins to work. 

The construction of such a quantum relativistic hadron theory has been given in papers 

P, m. 



5.9 The low energy theorems 



Let us check the existence of low energy theorems, in the bilocal approach ^ 



37| . The low energy theorems of hadron physics are independent of any model of hadrons. 



Consequently, they play a role of a pure tool to test the QCD-inspired models. These models 
(for example, the QCD sum rules [^, Nambu - Jona - Lasinio (NJL) model [^, QCD lattice 



calculations I^Tj, nonrelativistic quark potential models [|^. 11' EH) I'^'st^ the PCAC in hadron 
physics, induced by the spontaneous chiral symmetry breakdown, to the quark condensate. It 
should be noted that in all these quark models the definition of the local quark condensate is 
used. 



In the present section we discuss the validity of the notion of the local condensate in the 
quark potential model within the low energy theorems. Since a wide class of nonlocal potentials 
used in the model lead to equations that are unsolvable by an analytical method, we investigate 
the problem by means of the separable approximation. We calculate the amplitudes of leptonic 
and anomalous decays of mesons and consider the possibility of the generalization to heavy 
quark limit. 

As it was shown above the "constituent" quarks and their bound states are described by the 
following SD and BS coupled equations: 

f E(p)cos2??(p) =m° + i/g[cos2^?(p)] (la) 
1 E{p)sm2^{p)^{\p\) + li,\pqsm2^p)], (2a) 



MhUp) = Et{p)L,{p) - { [c(-)4-) + p ■ qs(-\s(-^] L,iq)} 
MhL,{p) = Et{p)L2{p) - i, { [c(+)c(+) +p- Uq)} . 



(299) 



Here mP is the "current" quark mass, Et — E1 + E2, Ei and E2 are the energies, respectively, 
of the quark and antiquark, and ^{p) is the solution to the SD equation, Li, L2 and Mh 
are the eigenfunctions and eigenvalue of the BS equation, respectively, identified by the wave 
functions and mass of pseudoscalar meson (H), and the following abbreviations are used: 

cW=cos[^i(j,)±^2(j,)], sW = sin[0i(j,)±02(p)l, j> q = ^, p = |p| , 



where ^(| p — q |) is the potential. We define the mass function as 

m{p) = E{p)cos2{p). (300) 
The solutions of the BS equation satisfy the normalization condition 

(301) 

In the chiral limit, these equations admit the solutions satisfying the Goldstone theorem. In 
other words, if a solution to the SD equation for the massless quark {mP — 0) exists, then the 
same solution satisfies the BS equation for a massless (M^ = 0) pseudoscalar meson too: 

L-(,) = L.(,) = ^^^f|M, (302) 

where is the pion leptonic decay constant having the experimental value 132 MeV. 

The latter equation is obtained from the amplitude of the decay tt — > in the chiral limit 
and normalization condition. In the potential model, the pion leptonic decay constant is defined 
as 

^'^ = ^1(0^'°'^^^^^^^^^^ 
where ^(p) = M^l±Mi_ is the pion mass. 



Notice that in the limit of heavy quarks << 1) the BS equation reduces to the 
Schrodinger equation 

2 

where 



i ' 

nil + m2 

whereas equation ( |303|) coincides with the nonrelativistic definition of the meson decay constant 

Thus, we can see that the quark potential model reflects the properties of the QCD at long 
and short distances. 

The separable approximation to the potential consists in the following factorization of the 
potential: 

nP-q) = |^/(|)-/(f) , /(0) = 1 , (304) 

where e is the coupling constant, / is the form factor, and L is the cut-off parameter. 

In this approximation one neglects the dependence of the spectra on the angular momentum. 
As a result, the SD equation takes the following simple form: 

m{p) = wP + mf[—) (305) 

or 



where the following abbreviations are used: 



m^m(O), <F>=j^,F ■ 



Let us consider the solution of this equation in the low energy limit, more exactly, when 
m°/L — ^ 0. In this limit, equation ( p06|) can be solved under the coupling constant 



r dx^££^=\ ~47rf r dxxfix)] , (307) 



where the approximation for the integral is valid as the dependence on 7 = m/ L is weak. 



We define the "averaged" quark Green function G{q) with the density function / as follows: 
« qq »= 4iV,ztr J J^^{f{q)G{q)) = -AN, < ^ >= -4iVc^m. (308) 
In the separable approximation, the BS equation also takes the form of the usual algebraic 



one 



< Ml 



m Aw? 



P(3)2 



1 - 



P(^) = |J<^> . (309) 



In the low energy limit, M^r << 4i?^(0), one obtains the following solutions to the equation: 

^^^^^ = T.mv^^'^^^)^^^) ' ^ ^ 

Now using eqs. ( pOlj ), (|306| ), and ( [308| ) we arrive at the following relation: 

- 4m° « qq »= F^M^ . (312) 

We see that this equation turns out into the well known low energy theorem of the local theory 
if one replaces the averaged Green function, << qq >>, with the quark condensate, < qq >. 
The quantities included in this relation can be estimated if the form of the potential and the 
current quark mass are given. 

To study the relation of the averaged Green function of the quark to the usual quark conden- 
sate, we consider the local NJL model type potential, f{x) = 6(1 — x), and use the conventional 
values for the free parameters, namely, 

L = 950MeV, m\ = AMeV, (313) 

then we obtain the following numerical estimations 

= 140MeV : = 132MeV ; 



(314) 



^ 7' = - 



«qq» = (-250MeV)3 ; = p = - , , 
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So, we see that these estimations are in agreement with the conventional values for the 
mass and decay constant of the pion as well as the quark condensate. We can conclude that 
equation ( P08| ) defines the quark condensate formed of the quark - antiquark pair due to the 
nonlocal interaction in the separable approximation. 

Let us now consider the mode of the pion decay, vr — » 77, induced by vector current. The 
amplitude of this process in the potential model is written as [|2T 



M{tt° — ^77) = < 77|VFgfp |7r° >= ^ ^ 



(27r)98PoA;ioA;2o 

(315) 



• 2iV,(e^(A:0e^(fc2)^(e^ - el)e,^^f,k^k§I{h, k. 



where e^{ki) is the photon polarization tensor (i = 1,2), are the quark charges, and Nc is 
the number of colors. The integral of this equation in the limit of zero photon momenta has 
the form 



I{k, = 0,k2 = 0) = J 



rfq m(q) 1 Ipr / n ^^(q) t i \^ 



(277)3 ^2(g)_M2 



(316) 



where the approximation of equal masses of the quarks is used, -E(g) = ^ + m{qY is the 
quark energy, P and fci, /c2 are the 4 - momenta of the initial and final particles, respectively. 
In the separable approximation the above integral takes the form 



J(A;i = 0, A;2 = 0) 



J 



J 



dx- 



(317) 



(318) 



1.02 



Vl [3(l + 72)3/2_ 

where for the numerical estimation the values of the parameters ( pi4| ) are used. 

So, the separable approximation for the integral J leads to the the Bell - Jackiw's theoretical 
estimation {J^-^ = 1) and the experimental result {J^'^'p ^ 1.04). 

Let us investigate whether one can self-consistently reproduce the low energy theorems for 
the pion decays using the conventional definition of the local quark condensate, namely. 



< qq>= ziVetr / J^,G{q) = ~2N^ J ^os^ P I)) 



(2' 



dq 

(2Tp 



(319) 



where G{q) is the quark Green function. It is easy to see that the low energy theorem can be 
obtained from equations ( |301|) , ( |302|) , and ( [319| ) if the "small" component of the bound wave 
function (L2(p)) in this energy scale has the behavior 



m 



R 



const . 



(320) 



However, this asymptotic function contradicts the Bell - Jackiw estimation. Indeed, substi- 
tution of 



into ([316|) leads to 

I{k = 0) 
Notice that by using 



Li{p = 0) ~ ^ , L2 



■J , J 







Li{q) 



dx 



X 



{l + x 



2^2 



4 



(321) 



(322) 



m 



instead of the value (|316| ) one would obtain ^7 = |. 

Thus, a straightforward application of the local quark condensate in the nonlocal potential 
model does not provide a self-consistent reproduction of the low energy theorems for the pion 
decays. 

The main result of the section is the check of the low energy theorems on the level of 
the effective bilocal meson theory and the generalization to heavy quark limit. This effective 
approach may be a reasonable for the detailed study of amplitudes and the Isgur-Wise function 
for heavy quarkonia decays (for example, B — > D*{p)lu) []5l| . 



5.10 U(l) - problem 

We have seen that the bilocal linearization of the four fermion interaction leads to an effective 
bilocal field t]m action 0, g in both QED ( |lTl] ) and QCD. 

This meson action includes Abelian anomalies in the pseudoscalar isosinglet channel [^, 
5U|, ^ (positronium r]M = VP' QED; and rjM = rjo-meson, in QCD). 

We have chosen the total-motion variable r]Mit) so that the effective action for the total 
motion of the the pseudoscalar bound state with anomaly term have the universal form for all 
gauge theories ((5-E'-D(3+i), QED^ij^i^, QCDi^^^i)) in terms of physical variables similar to ( |111|) 



where in 3-dimensional (5-E-D(3+i) the constant Cm is given by Eq. (|112|) , 



(323) 



C 



M 



Cp = — Svr 



2 I ILSch'- 
3/2 

nie 



in 1-dimensional Cm = 2v^; and in 3-dimensional QCD(^^^i^, 



C 



M 



X[A] is the "winding number" functional. In QED^^^i^ the "winding number" functional 
defined by Eq. ( |113|) describes the two 7 decay of a positronium. In QED^i^i) and QCD^-^j^i) 
these winding number functionals 



V/2 



Xqed(A^'^^) = 1^ / dxF^^e 



Nit) 



01 



v/2 
2 



2TtN 



3 I ^3 



(324) 



contain the independent topological variable. Recall that in QCD with the Wu-Yang monopole 
we obtained the normalizable zero-mode 



= NDf = NB^i^)- 



27r 



so that 



<iTT' 



J d^Df{<^)<^''^B^{<l>) = 1 . 



fin 



In QED(ij^i) and (5C-D{3+i) the effective action should be added by the topological dynamics 
of the zero mode with the actions 



where 



V/2 . 2 

-V/2 



l^QCD = \jdtj d'xGl = I dt^ 

V 



27r\2 1 

'QED - I — j y ^ 



'QCD 



27r\2 1 



.as J V<B'^> 

It is easy to show that the diagonalization of the total Lagrangian of the type of 

leads to additional mass of the pseudoscalar bound state in both QED{i^i) and (5CD(3+i) 

CI 



AM' 



IV 

In QED^ij^i) this formula describes the well-known Schwinger mass 

f<2 2 



AM' 



whereas in (5(71^(3+1) we obtain the additional mass of r]o meson 



Leff = llnl - v'omml + Aml)]V (325) 



Am/ = ^ = -4 ^ o • (326) 



'QCD' 

This result allows us to estimate the value of the vacuum magnetic field ( [232| ) in QCD 

< 5 >^ ' - 



2 



After calculation we can remove infrared regularization V — > oo. 

Thus, the Dirac constraint-shell formulation of gauge theories allows us to describe on 
equal footing the set of well-known results on anomalous interactions of pseudoscalar bound 
states in gauge theories including the anomalous decay of a positronium in QED(^r)=i+3), the 
Schwinger mass in QED(^e)=i+i), and U(l)-problem in QCD(^d=i+?,)- These results include also 
the zero probability of the two gluon decay of a pseudoscalar meson due to confinement as 
the destructive interference of the Gribov copies for gauge theories with the homotopy group 

T^{D-l){G) = Z. 

fil 



6 Conclusion 



Why does the non-Abelian gauge theory lead to the spontaneous chiral symmetry breaking, 
quark-hadron duality, rising potential, and color confinement? 

The formulation of the theory of strong interactions is completed, if the theory can explain 
and unify all working phenomenological schemes (chiral Lagrangians with the Goldstone mesons, 
the parton model based of quark-parton duality, the Schrodinger equation with rising potentials 
for heavy quarkonia, etc.) from the first principles (dynamic, quantum, and symmetric). 

The pure water of the theory of strong interactions can clean the muddy swamp of phe- 
nomenology if there are constructive explanations of the paradoxes: 

i) How to combine the observation of nonlocal gauge-invariant variables with the variational 
principles formulated for local fields? 

ii) How to combine the dependence of the Hamiltonian approach to quantization on the 
time axis with the relativistic covariance of S-matrix elements? 

iii) How to reconcile the presence of poles of color particle Green functions (that needs 
to ground the parton model and to obtain the Shrodinger equation) with nonobservability of 
physical states corresponding to these poles? 

The first two questions belong also to the consistent scheme of the description of the nonlocal 
bound-state sector (spectrum and interactions) in QED. Nevertheless, these questions become 
essentially actual in QCD where observables are only colorless bound states. In this case, the 
theorem of equivalence of different gauges |7^, |7^ on the mass-shell of elementary particles 
is not adequate to the physical situation of the description of observable nonlocal objects where 
local elementary particles are off mass-shell. For nonlocal bound states, even in QED, the 
dependence on the time-axis and gauge exists. All peculiarities of bound states (including time 
initial data, spatial boundary conditions, normalization of wave functions, time evolution) reflect 
the choice of their rest frame of reference distinguished by the axis of time chosen to lie along 
the total momentum of any bound state in order to obtain the relativistic - covariant dispersion 
law and invariant mass spectrum. 

Thus, the answers to the first two questions (i,ii) are the equivalent unconstrained system 
obtained by resolving the Gauss law constraint and the Markov- Yukawa prescription of the choice 
of the time axis parallel to the eigenvector of the bound-state total momentum operator. In 
this case, we can combine the relativistic and quantum principles of the dynamic description of 
gauge fields. 

The last question iii) belongs only to QCD. The answer is the topological degeneration of 
initial data in the non-Abelian theory that leads to the pure gauge Higgs effect of appearance of 
a physical vacuum in the form of the Wu-Yang monopole. This vacuum determines the effects 
of hadronization and confinement in QCD omitted by the conventional gauge-fixing method. 

We brought out that the physical vacuum leads to rising potentials, chiral Lagrangians, and 
an additional mass of the r^o-meson (as a consequence of the rearrangement of the perturbation 
series). The averaging over the Gribov copies of the topological degeneration leads to the color 
confinement and quark-hadron duality as a consequence of the destructive interference of the 
Gribov phase factors. This confinement is one of the consequences of the Dirac definition of 
measurable variables that keeps quantum principles. This confinement (omitted by all other 
methods of quantization) appears as pure quantum effect that gives the theoretical basis of the 
quark-hadron duality as the experimental method of measurement of the quark-gluon quantum 
numbers. 

Here we can recall the words by J. C. Maxwell in the Introduction of his A Treatise on Elec- 
tricity and Magnetism (Oxfrord, 1873): "The most important aspect of any phenomenon from 
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mathematical point of view is that of a measurable quantity. I shall therefore consider electrical 
phenomena chiefly with a view to their measurement, describing the methods of measurement, 
and defining the standards on which they depend." 
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